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ABSTRACT 



This thesis conducts an experimental investigation of 
the effects on separation of the closing angle of the stern 
section, the depth of submergence or proximity of the free 
surface, and the speed of a submerged body of revolution. 

The results of this investigation are then compared to the 
results obtained from two existing theoretical criteria for 
the prediction of the separation point on such bodies in 
turbulent flow. The point of separation is observed experi- 
mentally by photographing dye-flow patterns over the stern 
of the model at three speeds from .75 to 1.7 knots, at three 
depths of submergence corresponding to surfaced, near sur- 
face and fully submerged conditions for five different body 
stern shapes increasing in streamlining from a blunt hemis- 
phere to a finely shaped almost pointed section. Equivalent 
body profiles to the outside of the turbulent boundary layer 
displacement thickness are computed and potential flow pres- 
sure-velocity distributions based upon the equivalent pro- 
files are calculated and used in evaluating the separation 
criteria. 

It is concluded from the results that initiation of 
separation occurs at a nearly constant body angle on all 
bodies tested regardless of the degree of streamlining. Ex- 
perimental results show no apparent effect of depth of sub- 
mergence on the separation point and the effect of speed 
shows no clear cut pattern. The sensitivity of the variables 
used in the separation criterion on the predictions of the 
separation point are evaluated and recommendations are made 
for further studies and experimentation. 



Thesis Supervisor; Philip Mandel 

Title; Associate Professor of Naval Architecture 
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NOTATION 



Symbol 


Meaning Units 


C 


Non- dimensional coefficients 


E 


Error functions 


K 


Matrix functions of body 


K 1 


Matrix functions of prolate spheroid 


L 


Body length ft. or in. 


P 


One half the perimeter of the body profile ft. or in. 


R 


Gaussian weighting functions 


R s 


Reynold’s number, along body profile 


R x 


Reynold’s number, along x-axis 


S 


Surface area, over body ft.^or in. 2 

s 


U 


Total velocity along body ft. /sec. 


Uoo 


Free Stream velocity ft. /sec. 


V 


Velocity of body ft. /sec. or knots 


X 


Distance along x-axis ft. or in. 


Y 


Distance along y - axis ft. or in. 


ds 


Element of arc length along body 


g(x,t) 


Function of prolate spheroid whose ends coincide 
with body ends and intersects body at x - t 


g 


Acceleration due to gravity 


k 


Matrix functions of body 


k l 


Longitudinal virtual mass coefficient of given body 


k , (x) 


Longitudinal virtual mass coefficient of prolate 
spheroid at x - t 


n 


Reciprocal of exponent in boundary layer velocity 
profile 




Matrix functions of prolate spheroid 
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NOTATION (Cont ' d) 



Symbol 

P 

Q 

r 



r 



t 



u 

v 

w 

X 

y 

0 

0 

y 

A 

S 

8 * 

0 

X 

V 



Meaning Units 

Dynamic pressure lb/ft2 

o 

Stagnation pressure lb/ft 

Radius of body profile (calculation in. 

for corrected radius) 

Distance from location of point in. 

source on x-axis to point on body 
profile (calculatiQn of pressure 
distribution) 



Location of point source on x-axis, also 
equal to position of Gaussian abscissae on 
x-axis, normalized 

Component of velocity along x-axis ft /sec 

Component of velocity along y-axis ft/sec 

Axisymmetric potential function 

Distance along x-axis, normalized 

Distance along y-axis, normalized 

Angle between outward normal to body 
and x-axis (calculation of pressure 
distribution) 

Constant in separation criterion 

Angle between normal to body and y-axis or 
between tangent to body and x-axis 

Ratio of g*/g 

Boundary layer thickness 

Boundary layer displacement thickness 

Boundary layer momentum thickness 

Ratio of body length to body diameter 

Kinematic viscosity of water 



in. 

in. 

in. 
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NOTATION (Cont *d) 



Symbol 


Meaning Units 


f 


Normalized position of Gaussian abscissae 
based on body length of 2 and from -1 to 4 1 


P 


3 

Mass density of water slugs/ft 


9 


Velocity potential function 


9 

S 1 


Stokes stream function 

e 

Ratio of /£ in calculation of corrected 
body profile 


n 


Stream function of a unit source on the x-axis 


Subscript 


Meaning 


d 

f 

i 


Subscripts drag coefficient 
Subscripts coefficient of friction 
Gaussian abscissae index 


j 

k 

n 


Gaussian abscissae index 
Velocity is in knots. 

Refers to direction of body normal in derivation 
of pressure-velocity distribution 


n 


Index of iteration for solving equation of 
pressure-velocity distribution 


P 


Subscripts pressure coefficient 


Pi 


Quantity referred to flat plate 


r 


Subscripts residual resistance coefficient 


t 


Subscripts total drag coefficient 


X 


Quantity referred to x-axis 
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I - INTRODUCTION 



Separation of turbulent flow with its associated causes 
and effects has for the last half century been a subject of 
controversy in the field of naval architecture and ship de- 
sign. It is in this aspect of fluid flow that the naval 
architect has his greatest lack of knowledge. There is even 
considerable debate and controversy as to just what is known 
and what is not known on the subject of separation, its 
causes and effects. Little has been mentioned in the past 
regarding this aspect of fluid flow in references or texts 
on ship design. What has been stated has been conjecture 
based on limited theoretical calculations, sketchy experi- 
mentation, or recalled experience and observations. 

The causes and effects of separation and the para- 
meters relating to them are complicated and closely linked 
and must be fully understood before any definite design cri- 
teria of value can be established for application to ship 
and submarine forms. Separation of flow on such forms in- 
volves the three dimensional boundary layer and turbulence 
theories which have hardly been more than elaborated be- 
cause of their associated mathematical difficulties. Even 

! 

the mathematically simpler axial symmetric flow becomes in- 
volved with the problems of turbulent flow which is not fully 
understood or theoretically formulated for the relatively 
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simple case of two dimensional flow over a flat plate. 

The ultimate goal is to establish a valid, but simply 
applied, theoretical criterion to predict the separation 
point and a foolproof method of testing such a criterion 
with models to aid naval architects in improving ship de- 
sign. To do the latter the scaling effects involved in 
separation must be firmly established and understood as well 
as understanding the basic phenomenon of separation itself. 

The most important and difficult step in realizing the for- 
mulation of such a criterion is the determination of the sig- 

ta 

nificant parameters to be used and the quantl A tive effect of 
orch on the location of the point of separation. Faired surfaces 
with modest closing angles apperexitly retard separation but the 
limiting magnitude of the body closing angle or extent of 
fairing required to minimize separation is unknown. 

The effects of hydrostatic pressure and atmospheric 
pressure on separation in general, and in particular the 
effect on the limiting closing angle have been the subject of 
considerable debate. Saunders (j] # gives some information re- 
lating to these topics but this information has apparently 
never been checked or verified. The effects of the free sur- 
face on separation on a submerged body as it approaches and 
| ►/ 

penetrates the free surface is another area of considerable 

, / 

debate. Body speed and surface roughness are other possible 
parameters. 

* Numbers in brackets refer to works referenced in Appendix H. 



The object of this thesis is to investigate the ef- 
fects on separation of varying depth of submergence or proximi- 
ty to the free surface, body speed, and the closing angle of 

I 

the stern of the body. The results are necessarily qualitative 
in nature since the primary tool used in this investigation, 
dye-stream photography , is limited in the exactness to which 
the point of separation can be determined by noting the lo- 
cation of the familiar flow pattern of separation, i.e., 
abrupt departure of the streamlines of flow from the body 
surface. Originally it was hoped to relate the observed 
movement of the separation point, if any, with the independent 
variation of the parameters listed above to the form drag on 
the body. However, the results of the drag measurements on 
the model were unsatisfactory due to faulty or inadequate 
instrumentation and this could not be done. Two separation 
criteria , one developed by Stratford ( 2 ) and one described 
by Granville in ( 3 j) , are theoretically evaluated and the re- 
sults compared to the observed point of separation from the 
dye-stream photographs. This provides some insight as to 
validity of the criteria and the significance of the para- 
meters used therein. 

It is felt that the results of this thesis are a step, 
although incremental, toward determining what the significant 
parameters in separation of flow are and how they are related 
to the phenomenon of separation. 
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II - PROCEDURE 



A. General . 

The work presented here was performed in three phases. 

In the first phase model velocity and drag were measured 

experimentally and dye-flow photographs were taken of the 

/ 

model stern section for a systematic series of tests. In 
the second phase theoretical calculations were made to de- 
termine the point of separation on the body which could be 
related to that observed in the photographs. This phase in- 
volved the calculation of the turbulent boundary layer dis- 
placement and momentum thickness, the pressure-velocity dis- 
tribution over the resulting corrected body and the use of 
these quantities in two separation criteria for turbulent 
flow. The third phase was the analysis of the results of 
phases one and two to determine and relate, if possible, the 
effects of the independent variation of speed, depth of sub- 
mergence, and stern shape on the observed separation point 
location on the body as well as to compare and relate the ob- 
served point to the calculated point of separation. 

B. Experimental Phase. 

In this phase experiments to visually determine the 
point of separation on the model and to measure the total 
drag on the model were conducted at the Massachusetts Institute 
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of Technology Ship Model Towing Tank under the supervision of 
the Department of Naval Architecture. Drag and velocity data 
and dye-flow photographs were obtained for the model at three 
speeds, at three depths of submergence corresponding to sur- 
faced, near surfaced and submerged conditions, and with five 
varying shaped tail sections which modified the closing angle 
of the model's stern. 

General Description of Equipment. 

\ 

Figure I is a sketch of the model showing the variation 
in stern section and the waterlines of the three depths used. 
Figure II shows a view of the towing tank and towing carriage, 
while Figure III shows the carriage control station and in- 
strumentation. 

The model, #4893, was constructed at the David Taylor 
Model Basin, Car derock, Maryland, and is & blunt nosed circular 
cylinder consisting of three main parts: a blunt solid ma- 
hogany nose section, a long hollow aluminum parallel middle- 
body, which contains the dye system, and five solid mahogany 

stern sections of varying shape and length. To insure opera- 

trip 

tion in turbulent flow a 1/32 inch diameter/wire was placed 
around the nose of the model 2.5 inches from the forward end. 
The offsets and other characteristics for the model and stern 
sections are given in Appendix A. Figure IV shows the model 
components . 
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FIGURE I . MODEL AND STERN SHAPES 




Figure II - MIT Ship Model Towing Tank 




Figure III - Left to right: carriage control console, 
digital voltmeter and speed measuring equipment, Sandborn 
recorder. 
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Figure IV - Model and five stern shapes 




Figure V - Internal view of model parallel middlebody showing 
support plate, dye tubes, dye manifold and brass dye ring. 



8 




The dye system emits dye from 1/16 inch diameter holes 
spaced every 1/4 inch apart around a circular ring located at 
the after end of the parallel middlebody. This system is shown 
in Figure V. Dye was fed to the model from a tank mounted on 
the towing carriage through 1/4 inch plastic tubing. A solu- 
tion of water and potassium permanganate was used for dye. 

Each of the five stern sections can be attached to l^he mo- 
del immediately aft of the dye ring. When the model was in 
motion dye released from the clye ring would stream over the 
stern section permitting observation of the flow. 

The model was rigidly mounted to the towing carriage 
through a bracket which allowed operation at any depth within 
the limits of the 18 inch towing street. Figure VI shows 
the model mounted on the carriage. The towing carriage rides 
on an overhead rail and is driven by a variable speed elec- 
tric motor driving a soft rubber wheel in contact with a 
similar side rail. 

i. 

Figure VI also shows the drag measuring dynamometer 
located between the 4:1 ogive towing strut, connected to the 
model internally, and the 2^ inch diameter hollow-cylinder 
which was connected to the mounting bracket. The dynamometer 
measures the axial force between the carriage and the model 
by means of a differential transformer. 

A length of 32 feet was used for each run. The length 

of run was limited to this length due to a portion of the 

1 
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Figure VI - Model mounted on carriage shov/ing towing dynamometer, 
support bracket, dye tubes and towing strut. 




Figure VII - 4x5 view camera and strobe light arrangement 
at tank window. 
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carriage rail beyond this length being badly distorted and 
thereby causing undue carriage vibration. The length of run 
started at the beach end of the tank and extended to the end 
of the first of the 10 foot by 4 foot glass windows located 
in the tank's side wall. A 4x5 view camera was set up at 
this window along with two high intensity strobe lights which 
were tripped by the carriage striking a contact on the over- 
head rail as it reached the proper position in the window near 
the end of the run. A white } light reflecting^ background was 
installed behind the model because the distance to the far 
side of the tank was too great to achieve satisfactory expo- 
sures. Figure VII shows the photographic set-up that was 
used. Each stern section was marked at 1/2 inch intervals 
along its length as an aid in determining the position of 
separation from the photographs. 

The carriage speed was set by the manual speed control 
dial on the carriage console since the automatic speed con- 
trol was inoperative. Speed data for each run was recorded 
on tape with a Sandbarn Type 150 recorder and type 1300 preampli- 
fier, Drag force was read out on a type 405 Hewlett Packard 
D.C. digital voltmeter with input from the type 1100 Sandborn 
preamplifier. The voltmeter's purpose was to attenuate the 
high frequency components of the carriage vibrations sensed 
by the dynamometer as drag fluctuations and to allow a volt- 
age proportional only to model drag to be read out. A more 
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detailed description of the towing tank and associated equip- 
ment, model dye system and dynamometer is given in Appendix B. 

Run Procedure 

Two runs were conducted for each condition of stern 
section^ depth and speed. On each run a photograph was taken 
of the dye-flow over the stern section while speed and drag 
were recorded simultaneously on the Sandborn recorder and 
digital voltmeter. The tank was allowed to settle between 
each run. Polaroid pictures were taken occasionally to check 
exposure and dye-flow definition Only five of the dye holes 
in the dye ring were used, two each on the top and the bottom 
and one on the axis of the body. In turbulent flow the dye 
present from the other holes tended to obscure the separation 
point. For laminar flow the use of the entire ring was satis- 
factory. Figure VIII shows both turbulent and laminar flow. 
The unused holes were masked off. All screw holes, etc, on 
the model were filled smooth with a waterproof putty to insure 
as little flow disturbance as possible. 

The model speed was easily obtained from the Sanborn 
trace and a speed calibration curve; however, the drag could 
be obtained from the digital voltmeter only by averaging the 
reading over the run. This was done by recording about ten 
values of the voltmeter reading during the last third of the 
run and averaging. This average reading was then compared 
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to the "zero" drag reading taken prior to each run to determine 
a voltage difference reading proportional to drag. A calibra- 
tion curve was determined by applying a known static load to 

the model and observing the voltmeter reading. Using this curve 

** 

and the voltage difference reading drag force was obtained. 

C. Theoretical Phase . 

In this phase a theoretical position of separation on 
the body was established after correcting the body for the ef- 
fects of boundary layer growth. Throughout this phase it was 
assumed that the flow was turbulent from the nose of the body 
and that the 1/7 power law holds for the velocity distribution 
within the turbulent boundary layers. 

Calculation of Turbulent Boundary Layer Characteristics . 

As indicated in the general discussion of the procedure 
the value of the displacement thickness on the body and mo- 
mentum thickness on the body are necessary. Eckert [4] has 
evaluated the effects of transverse curvature on the turbu- 
lent boundary layer on a circular cylinder in compressible 
flow. These effects are the same as reported by Yu ;(5) 
and Landweber (jf) , which are: 



** The drag measurements were unsatisfactory. Appendix E 
has a more detailed discussion relating to these measurements. 
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(1) Boundary layer thickness and displacement thickness 
are less than on the flat plate, 

(2) Momentum thickness and local and mean skin friction 
coefficients are higher than flat plate values. 

(3) The effects of (1) occur when the ratio of o/p^O.l. 

(4) The friction coefficients deviate from flat plate 

values only when ratio 1,0. 

Eckert develops simple equations which relate the value 
of , and © on the body to values of the same parameters 

on a flat plate. The assumption is made that the velocity 

profile within the boundary layer conforms to the l/7th power 

/ \ 1/7 

law distribution of U = fY] and that Blasius 1 Law relating 

Uoo \<5/ 

wall shear stress as inversly proportional to the fourth root 
of the boundary layer thickness are uneffected by the curva- 
ture. The assumption is also made that the boundary layer is 
turbulent from the beginning and it starts from the nose of 
the body with zero thickness. 

The relationships developed for incompressible flow are 
given below and the reader is referred to [4] for the deriva- 
tion details which are straightforward. 

Relationship for boundary layer thickness 

6 - A _\ 4/5 ( i ) 

T ~ pi \ 1+ 4/3:7 

where } 

£ Pi = 0.38 X (ia) 

(Rx) U - 2 

Although (1) is implicit in S> it can be easily solved by using 
£ pi as a first try and iterating the results. 
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Relationships for Displacement thickness. 




§_ 

r 




4/s 



where , 

A^pi 



6 p i 
Tpi 



i_ 

77+1 



V 1 +Ti 



(2b) 




An ( ZL 2 

h(2t>+1) 



1 

2 7 ? 4 1 



(2c) 



IX - T'A' orv 
lloo ’A ' 



Relationships for Momentum thickness. 




4 



£ 

2r 



Ao 



n-I 



where , 

O ' '' - 7b To 

7b +2 



7^ 9 pi 

-7L+1 (rt+2) TpI 



(3) 

(3a) 






2n.pi 




( 77 , 4 - 2 ) 

1*72+1) 



72/ (3b) 

( 77 41) ( 2 tv4 1 ) 



The value of ,v > * calculated from (2) was applied as a 
correction to the body surface after the manner described in 
(73 as is shown in Figure X taken from (Y) . 
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_ j j xS3>oy 

Figure X - Equivalent Body 
Profile with Displacement Thickness 

Thus the corrected body will have a radius given by 
V corr * V + Cos X £ (4) 

For use in the calculation of the pressure velocity dis- 
tribution discussed in the next section this radius was norma- 
lized as follows: 

3 ° corr - l^corr (4a) 

norm body length72 

An IBM 7090 computer program, Program II, was written to 
evaluate equations (1,2, 3, 4 & 4a) at 16 points on each of the 5 
bodies for 9 speeds. These 16 points correspond to the loca- 
tion of the 16 ordinates required for the Gaussian 16 point 
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quadrature rule use in the determination of the pressure veloci- 
ty profile. The values of body radius and the angle 6 were 
picked off at each of the 16 points from full scale drawings of 
the 5 bodies prepared by the authors. A tabulation of this 
data is contained in Appendix A. Details of the computer pro- 
gram are in Appendix C, 

Calculation of Pressure Velocity Distribution 

Potential Flow Pressure and Velocity Distribution was 
calculated according to a method developed by Landweber [8j . 

Two types of solutions are given for obtaining the steady, ir- 
rotational, axisymmetric flow of an inviscid, incompressible 
fluid about a body of revolution. The first method involves an 
iterative type solution for Fredholm integral equation of the 
first kind and is not used in this work. The second method ob- 
tains the velocity directly as the exact solution of an inte- 
gral equation of the first kind by an arithmetic procedure. 
Using the second method as outlined in QQ the basic problem 
is to solve for the velocity potential 0 which is dependent on 
x, y and satisfies Laplace’s equation in cylindrical coordinates 




where the equation of the normalized body profile is in 
the form: y^ - -f (Figure XI shows method of normalizing the 

body) 

and u “ - o>0 4/ - - c)0 

<57 “57 
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> 7 - 



Figure XI - Normalized Body Profile 
There also exists a Stokes stream function ^(x,y) 
which satisfied the equation: 

c> 2 ^ + 1 bP ( 6 ) 

dx 2 dy 2 y &y~ 

For uniform flow of velocity 11 parallel to the x-axis 
0 “ -£Zx and^- - hUy^ 1 at infinity 
The following boundary condition on the body must also be 
satisfied: 

/ Q0 \ = - V cos/g 
KcPto )s 

Where V - velocity of the body and /O- the angle between the 
outward normal to the body and the x-axis. 

Derived from Green's Theory the equation: 

Jf 0 &“ dS " ff aS- dJf (7) 



19 



where 0 and w are any functions harmonic in the region exterior 
to a given body and vanishing at infinity, is used to obtain 
an integral equation for axisymmetric flow about a body of 
revolution. 

Substitution of: dS - 2 ~TT y ds and 

d0 - - sin # into (7) gives: 
dn 



70 ds 



y w sin # ds 



( 8 ) 



y - ordinate of meridian section of the body which is 

moving in the negative x - direction with unit velocity, 
ds- element of arc length along body 
X - angle of the tangent to the body with x - axis 
2P- perimeter of a meridian section 
With \a I as an axisymmetric potential function and 
(x,y) the corresponding stream function 



y dw 
dn 



ds 



and 



P 

y0 dw ds 
dn 





ds 



If UL is total velocity along a body when a stream of 

unit velocity is superposed in the positive x - direction 

H -~d0 + cos# 
ds 

Since dx - ds cos# and dy - ds sin# equation (8) can 
be written as: 
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/P 

£/^ds=j ( dx - yw dy) - 

Jr 



p 

o 



(9) 



but y ■ 9w - - since w and are corresponding potential 

57*- Sy 

and stream functions and dx - ywcly is an exact differen- 
tial defining a function S~2 (x,y) such that 



0-n - tp 



d£2 - - y w 

< 5 y 



From this and y bw “ we obtain 

dy Vi- 

' c) + 3 2 X 2 + i^ 9X2 

y Ty 



a 



3 y 2 



which is identical to (6) and is the equation satisfied by 
Stokes stream function. It can also be verified that w and ^ 
are corresponding axisymmetric potential and stream functions 



which satisfy the equation: 9^-' - - y 

a* a y 

Writing (9) in terms of _CL 




QX2 ds ” ^X2 ~ 0 < 9-^2 \ 



a* 



/ 



9^- y 90 

9y c?x 



( 10 ) 



o 



If /I is the stream function of a point source of unit 
strength situated at x - t cm the x - axis 



Si 



i 



- 1 + x-t 

Jl 



Jl - jj(x-t) 2 + y 2 ] (11) 



Then - 2 

d x a 



y , and if y vanishes at both limits, 
n3 

P 



( vSl- 0 dSl ^ 
v dy^ t 



- 2 
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1 



( 12 ) 



(11) now becomes: 



L 




J o 

which is an integral equation of the first kind in which the 



where q is the stagnation pressure. 

Landweber uses an iteration procedure to solve ( 12 ) 
as follows: 

If the iterations are expressed in terms of error functions 
E n (t) at the point x = t , first obtain 



unknown function is ££ (x) and the kernel is y^^jp' ^ en u 
is determined the pressure distribution is given from 
Bernoulli’s equation. 



P_ - 1 - (u 2 + x/2) = 1-U 7- 



q 



E x (t) from E n 




by using the first approximation IJ. (t) - (l+k) cos X (x) 
where k., = longitudinal virtual mass coefficient. Then sue 



cessive iteration £3, E3 



E are obtained from 



n 



E 




dx 



whereto ♦ 0^ are nose and tail abscissaeof the body , and 



finally U n +j(t) is obtained from 



n 



n 



Un+i (t) " (t) -+ cos ^ (t) £1 E ± (t) ( 13 ) 



i = 1 
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To evaluate integrals of the form 



E (x) y (x) 

/t 3 

use k (x,t) - -f (x) ( 14 ) 

-t) 2 4 f (x) J^2 

k' <^t) - g (>, t) (15) 

[( *-t) 2 + g()Oj72 

where y 2 *» f (y ) is the equation of the given body profile 
and y 2 - g(x^ t)) is the equation of the prolate spheroid 
whose ends coincide with the ends of the given body, and 

f 

which intersects the body at x « t, i.e #) 

g(*, t) «/(t) ,(x-y 0 )(xi -x) (16) 

(t- x 0 )(Xi - t) 

I . 

For the spheroid the length - diameter ratio A is given 

by 7 l 2 - (t - xp) (x, - t) (17) 

/ (t) 

and from this the virtual mass coefficient (t) can be ob- 
tained from: 



/Unpi+y^r^i) - ta -1 

7-2 Ji - >ln(/W 2 _i7 



(18) 



E i (t)is then obtained from 



Ej, (t) 



and finally 



1 - 1+k, 



[k (xvt) -K 7 (x.t)^ dx - 1 + ki 

L U (t ) 



(19) 



E n+1 " ? (t)+ki E n (t) I k(X, t) rE n (x)-E n (t)1dx 
I + k / L J 



'Vo 



( 20 ) 
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In using Landweber's method the velocity is obtained 
at 16 points along a body, from t ■ -1 to t ■ + 1, correspond- 
ing to the Gaussian values ( ^ ) for the 16-point guadrature 
rule. In applying the rule the integrands in (19) and (20) 
are evaluated at the 16 Gaussian abscissae x j = for each 
of the 16 values of ti«* f ± . Gaussian weighting functions 
are used to obtain K j i = R j k (Xj, tj.) and K j i - Rjk 
(Xj } tj^) which are used in evaluating (19) and (20).U n (t) is 
then given by (13) and ‘p/^. by Bernoulli’s equation in the form 
P/q " 1 -TJ n 

In the application of Landweber’s method of determining 
the pressure-velocity distribution on a body of revolution to 
the models used in this thesis the corrected body profile as 
calculated in the previous section was used. The function 
y2 -^(x) for the normalized corrected body profile was not 
determined because even if it is determined the end result is 
still a point-wise distribution of^J/Lk>. Instead of determin- 
ing the polynomial describing the body profile the values of 
the normalized corrected body radius at each Gaussian ordinant 
computed in Program II were used. The angle, # , was again de- 
termined from full scale drawings of the model and stern sections. 

An IBM 7090 computer program, Program III, was written to 

i 

determine the pressure-velocity distribution on the body cor- 
rected for boundary layer effects using Landweber’s method. The 
program calculates Cp, VL/\Joo» absolute velocity and gage pres- 
sure on the body at 16 points on each of 5 bodies for 3 different 
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speeds at each of 3 different depths of submergence. The 9 
speeds used are the same as those used in Program II and are 
the average of the speeds listed for each run and its dupli- 
cate, i.e., Run 1 and Run 1A, in the tabulation of experimen- 
tal data in Appendix D. Thus the outputs of the two programs 
are a compatible package for application in the separation 
criterion to be evaluated. This output is in turn compati- 
ble with the dye-flow photographs since the same conditions 
apply. A test solution for Program III was run on the ex- 
ample worked out in [ 7 ~] . 

Separation Criteria . 

Stratford [_ 2 ] has developed a rapid method for the pre- 
diction of flow separation in turbulent flow from an approxi- 
mate solution of the equations of motion. The equations are 
integrated by a modified "inner and outer layer solution" 
technique developed by Stratford in 1954 for laminar flow 
boundary layers. The final solution is given by equations 
(21) through (24) which relate X, Cp, and dCp/dX at the sepa- 
ration position and may be applied to any given pressure dis- 
tribution. The algebra and procedure for obtaining the sepa- 
ration condition are clearly detailed in [2] and only the re- 
sults will be given here. 

The resulting implicit equation for the separation po- 



25 



sition, X, is 



i 



(2Cp)^ n “ 2) (x - 1 . 06 (l(T 6 Rs) 



- 1 

TO (21) 



for Cp 



n+1 



Where: R g is the Reynolds number based on value of distance 



X and freestream velocity, Uoo. 

n = log 10 R s 

and - 0.66 for d 2 p <. 0 

dX 2 

ft - 0.73 for d^p > 0 

dX 2 



( 22 ) 

(23) 

(24) 



with the resultant value of d p being that immediately prior 

dJT 



to separation. This formulation has a range of uncertainty 
of 10% due to its inability to fully account for d 2 p/dX in 
the case of turbulent boundary layers. However, the rapidity 
of the method makes this range of uncertainly acceptable in 
this case. Granville in [3] has brought forward another even 
simpler separation criterion It was apparently originally 
developed by Bu*i in 1931 at Zurich. The criterion states 
that separation will occur if the following relation holds: 

© Uoo du L < -0.006 (25) 

u * L dX Uoo ~ 

Where; 9 = momentum thickness on the body at separation 

point (Calculated using equation (3^ ) and L - body length. 

) 

since L and Uoo are constants for each body at a parti- 



cular speed 



U 



( 26 ) 



Equation (25) can thus be written 




Unfortunately the information regarding this criterion is 
limited to that obtained by the writers from Course 13.012 in 
the Department of Naval Architecture at M.I.T. and the meager 
yet reliable reference to it in [3] . A range of uncer- 
tainty for the criterion is not known. It is intended to be 
used as a criterion to be used in scaling separation from 

model to ship scale. 

( 

With the velocity profiles, pressure coefficients, and 
momentum thickness values on the body determined as outlined 
in the previous sections it is possible to evaluate both of 
these criteria and see how they will compare with the experi- 
mentally observed points of separation. This was done by 
hand calculation for a few cases. Sample calculations are given 
in Appendix F. The results of such calcuations are in the 
RESULTS section of this thesis. 

D. Analysis Phase. 

The data obtained from the experimental phase relating to 
drag was analyzed by utilizing an IBM 7090 computer program, 
Program I, to reduce the data to coefficient form. The details 
of this computer program are given in Appendix C. The analysis 
of the experimental drag data contained in Table D-I of Appen- 
dix D showed the data to be unreliable and faulty. Appendix E 
contains some of the results obtained from the computer program 
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and a discussion relating to these results. The second step 
in analyzing the experimental data was to study the dye-flow 
photographs and determine from them the point of separation 
and the body angles with respect to the body axis at the point 
of separation. The point where the dye-flow appeared to a- 
bruptly depart from the body profile was taken to be that of 
the point of separation. Only the flow along the keel of the 
model was used in determining the location of the point of 
separation because it was felt that the presence of the towing 
strut forward would influence the flow over the top of the 
model. The authors each made independent determinations of 
the separation position on each run. Differences which arose 
in interpreting the location of the separation point were 
then compromised. The determinations were also made prior to 
the calculation of the separation position from the two sepa- 
ration criteria so that the location of these theoretical 
points would not influence the interpretation of the photographs. 
In this manner it was felt that as high as possible degree of 
precision would be obtained in interpreting the photographs. 

There was, as was to be expected, considerable judgment in- 
volved in establishing the separation point in several of the 
pictures. A sampling of typical photographs from which the 
separation point was determined is shown in Figure XII. The 
location of the observed separation 'point for each run is con- 
tained in Table D2 in Appendix D. The angle between the body 
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Run U 
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Run 7 

Ac 




Run 13A 
Body 1 




Run 20A 
Body 2 



B 

Figure 





Run 3 




Run 5 A 




Run 6 




Sequence of Runs on Body 3» 
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Run 12A 
Body 1 




Run 35A 
Body U 







Run 29 
Body 4 




Run 38A 
Body 5 



. Runs on Other Bodies. 

XII - Typical Dye-flow Photographs 
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surface and the body axis at the point of separation was de- 
termined graphically by plotting the separation point loca- 
tion on full scale drawings of the bodies. The movements of 
the point of separation were analyzed to see what, if any, 
correlation could be made between these movements and changes 
of speed, depth of submergence or body shape. The angle 
at which separation occurred .was similarly analyzed. 

The third step in the analysis was to calculate the 
point of separation by each of the two separation criteria 
using the results of computer programs I and II. Samples 
of the output from these programs are contained in Appen- 
dix G. Sample calculations for each of the criteria are 
given in Appendix F. During these calculations the criti- 
cality of the various parameters in the criteria were 
noted. The results of these calculations were compared 
to the observed point of separation noting location and 
angle of body surface at the point of predicted separation. 
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HI RESULTS 



The results of the procedures described in the pre- 
vious sections are presented in this section. These results 
are discussed and conclusions are drawn from them in the 
next section of the thesis. 

The results are presented primarily in graphical 
form. Figures XIII through XVII show the relationship be- 
tween the pressure-velocity distribution over the body and the 
location of the observed separation points. Figures XVIII 
through XXII show the location of the experimental and calcu- 
lated points of separation. The body profiles corrected for 
boundary layer effects are also shown here. Table I con- 
tains a breakdown of the experimentally observed separation 
points on each body by depth of submergence, and speed. The 
body angle at separation is also shown. 
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TABLE I 



Results of Studies of Dye Flow Photographs and Separation Criteria 





BODY # 1 


X/ 








Code 




Angle 


/ L 


Points (Speed-Depth) 






A 


- 


240 




5 @ LM, 2 3 LD, 2 @ MD, 2 




LS 


B 


- 


32° 


.970 


3 @ MM, 2 i HM, 1 e MS 






C 


— 


420 


.978 


2 e HD, 2 g HS, 2 e MS 






E 


— 


60 


. 948 


BURI CRITERION 






F 


- 


47° 


.982 


STRATFORD CRITERION 








BODY #2 


X / L 








Code 




Angle 


Points (Speed-Depth) 






A 


- 


“25 s ” 


7557 


1 S Ld, 1 @ LM 






B 


— 


30° 


.976 


2 ® MM, 1 € HM, 1 fe LD, : 


2 @ LS 


C 


— 


36° 


.983 


2 e MD, 1 4 HD, 1 e MS, 1 




LM 


D 




42° 


.989 


2 @ HS, 1 @ MS 






E 


— 


13° 


.930 


BURI CRITERION 






F 


— 


38° 


.986 


STRATFORD CRITERION 








BODY #3 


X / 








Code 




Angle 


x /l 


Points (Speed - Depth) 






A 


- 


23 5 


.957 


2 4 MS, 2 € LS, 1 e HS 






B 


— 


28 o 


.975 


2 4 LM, 2 € MM, 1 e HM, 1 


@ 


HS 


C 


— 


33 


.982 


2 4 HD, 2 4 MD, 2 4 LP 1 


e 


HM 


E 


- 


10° 


.905 


BURI CRITERION 






F 


— 


46° 


.992 


STRATFORD CRITERION 








BODY #4 










Code 




Angle 


x /l 


Points (Speed - Depth) 






A 


— 


26 ° 


.983 


2 4 LD, 1 @ MD 






B 


— 


38° 


.990 


2 « LS, 2 4 MS, 1 @ HD, 1 


4 


MD, 


C 


— 


49° 


.996 


2 4 HM, 2 ^ MM, 2 ^ HS, 1 




LM 


E 


— 


6° 


.881 


BURI CRITERION 






F 


— 


36° 


.989 


STRATFORD CRITERION 








BODY # 5 


X , 








Code 




Angle 


/l 


Points (Speed-Depth) 






A 


- 


23“ 


.978 


All speeds at deep depth 






B 


- 


30° 


.988 


All speeds at medium depth 




C 


- 


60° 


.998 


All speeds at surface 






E 


— 


5° 


.831 


BURI CRITERION 






F 


- 


23° 


.979 


STRATFORD CRITERION 







Separation points are plotted on body profiles in 
figures XVIII through XXII. 

Point designations are the same as outlined in Table D-II 
of Appendix D; e.g. 2 @ HM means two points were plotted at 
that location and angle on the body corresponding to two runs 
at the high speed - medium depth condition. 
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FIGURF XT7T1I SEPARATION POINTS . BODY 
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IV DISCUSSION OF RESULTS 



Accuracy of Observations 

The accuracy achieved in determining the location of 
the separation points from the dye-flow photographs was 
limited by the ability of the authors to pick off the point 
where the flow abruptly departed from the body surface and 
in their ability to interpolate between the £ inch markings 
on the bodies. In some cases the point of separation was 
partially obscured by excessive amounts of dye while in 
other cases, especially at the higher speeds, a portion of 
the dye appeared to be carried beyond the separation point 
by flow outside the boundary layer. In the surface condition 
there was some evidence of dye being drawn up into the area 
aft of the body toward the free surface. The accuracy in 
observing the point of separation in turbulent flow can not 
be equal to that obtainable for laminar flow, but improvement 
in this accuracy could be gained by using a larger model and/or 
using a more sophisticated means of photographing the flow 
such as telephoto lens or underwater motion pictures to ob- 
tain greater definition of the flow about the body. 

Effect of Body Closing Angle 

In examining the positions of the observed points of 
separation on the various bodies as shown in Figures XVIII 
through XXII it was noted that the points occurred in groups 
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located at three or four positions on the bodies. Table I 
shows the details for these groupings. The center of each 
group was used as the measuring point of the body angles. The 
degree of accuracy obtained in measuring these angles is a 
function of the accuracy achieved in the initial location of 
the point of separation already discussed and the errors en- 
voi ved in measuring the angle from the full scale drawings. 

The two factors result in a possible error of - 3 degrees. 

i 

The initial group of separation points on all five 
bodies tested occurs at a body angle of between 23 and 28 

I 

degrees, averaging about 25 degrees. There is no apparent 
orderly correlation among the five bodies to account for this 
variation except that it is due to experimental error. This 
results in the conclusion that the initial point of separa- 

I 

tion occurs at approximately the same angle for all bodies 
regardless of the amount of streamlining. A critical angle 
of 25 degrees thus appears to be established for initiation of 

I 

separation. The rate at which the body angle is changing does 
not appear to effect this critical angle because bodies 1 and 5 
which have vastly different rates of change of body angle have 
nearly the same critical angle. A systematic check over a 
greater range of body shapes, speeds and depths of submergence 
is warranted to check this conclusion. 

The maximum value of the angles in the region over 
which the separation points occurred varied from about 42 degrees 
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to 60 degrees showing that the observed point of separation 
did fluctuate. At these higher angles, however, the speeds 
were generally toward the high end of the speed range and the 
fluctuations may be caused by dye being carried aft by flow 
outside the boundary layer, especially on body 5 where the 
point of separation was most difficult to determine. Be- 
cause of this possible effect these angles are not con- 
sidered as significant as the initial separation angles. 

The above comment, notwithstanding, these higher angles do 
point up that under certain circumstances separation of 
flow may not occur until the critical angle has been ex- 
ceeded by a considerable amount. 

Evaluation of Separation Criteria 

Buri's theoretical separation criterion equation (26), 
consistently predicts the separation point well forward of 
both the experimentally observed points and the point of 
separation established by Stratford’s criterion, equation 
(21). For Buri's criterion the body angle values all fall 
between 5 and 13 degrees with no apparent consistency of 
variation, except that when body 1 is eliminated the angles 
decrease from body 2 to body 5, thus indicating that the 
separation point moves forward with increase in streamlining. 

* i i 

This is inconsistent with the experimental results and the 
general concepts relating to separation which maintain that 
streamlining should decrease the extent of the separation zone, 
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not increase it. This criterionis very sensitive to slope of 
the curve of ^Uoo versus X ^L as plotted in Figure F-I, Appen- 
dix F, and relatively insensitive to model speed that enters 
through 0 in (26). The criterion predicts separation to occur 
lj.ust as the value d^Uoo)^/d ( X// l) becomes negative. A slight 
error in plotting this curve or in using it to solve (26) 
could result in a larger error in the separation point; how- 
ever, the magnitude of the possible error, -0.2 inch, is not 
sufficient to improve the agreement between Buri and the ob- 
served points. 

On the other hand, Stratford’s criterion gives much 
closer agreement with the experimentally observed points of 
separation, as it consistently predicts separation either with- 
in or slightly aft of the observed region of separation. The 
values of the body angles at the predicted separation points 
vary from 23 to 47 degrees. Again there is no apparent con- 
sistency of this variation. Depending on the shape of the curve 
of Cp versus X, as plotted in Figure F-II, Appendix F, this 
criterionis sensitive both to the value of Cp and the value 
of [ p’/dx, and also insensitive to model speed that enters 
through R and n in (21). Where the slope is steep the value 
of Cp is the controlling factor while for a shallower slope 
both Cp and d(CpVdx control. The criterion essentially pre- 
dicts separation to occur very shortly after Cp becomes posi- 
tive. Again errors in plotting or interpreting Figure F-II 
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would not effect the predicted point of separation by more than 



-.0.1 inch. 

It can be concluded that Buri's criterion did not give a 
satisfactory comparison to the observed point of separation 
while Stratford's agreed much better although it generally 
indicated that the body angle at separation was of a larger 
value than that observed. It is to be kept in mind, however, 
that these criteria were analyzed on the basis of the calcu- 
lated profiles in Figures XIII through XVII which were de- 
rived for modified potential flow. It is therefore recom- 
mended that an experimental determination of velocity-pres- 
sure profiles be made on the bodies and that these be applied 
to the separation criteria before a final judgment is made. 

Effect of Proximity of Free Surface 

Table I yields no apparent consistent variation of the 

experimental points of separation on the various bodies with 
a change in operating depth. 

Computer Program III calculates gage pressure or dy- 
namic pressure plus hydrostatic pressure due to depth of sub- 
mergence on each body. The variation in pressure at the keel 
is from about 19 lbs. per square foot at the surfaced condi- 
tion to about 93 lbs. per square foot at the deep depth which 
should i c sufficient to show the effect of change in depth on 
the observed location of the separation point if in fact such 
an effect on separation does exist. Both of the separation 
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criteria are completely independent of depth since it does not 
enter into any of the variables contained in equations (21) or 
(26). It is concluded from the experimental results that hydro- 
static pressure over the range investigated has no recognized 
consistent effect on the location of the separation point. 

It is recommended that further experimentation be 
conducted with larger variations in depth. Perhaps the use 
of a water tunnel or propeller tunnel would be better than 
a towing tank for this type of experiments since a large 
static pressure could be applied to the model in the test 
section. More typical ship shapes rather than body of re- 
volution forms should be used in conducting dye-flow experi- 
ments in the air-water interface to show the effects of 
draft on separation over a given form. 

Effect of Speed 

Bodies I through 4 show values of the body angle at 

the separation point between 23 and 28 degrees for the low 
speed range, between 30 and 37 degrees for the medium speeds 
and 33 to 60 degrees for the high speeds. Because of the dif- 
ficulty in locating the separation point for body 5, no rela- 
tion between speed and separation point could be seen for that 
body. The results for bodies 1 through 4, however, show that 
the separation point appears to move aft on the body as speed 
increases. The reason for this effect could be due to the 
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fact that dye was being carried back along the body outside 
the boundary layer and misinterpretation of the photographs 
resulted. At the speeds used the Reynolds' Numbers were from 

5 

3.5 to 7 x 10 . This is within the doubtful transition region 
between laminar and turbulent flow and the possibility exists 
that the lower speeds were at the laminar end of the transi- 
tion, although an examination of the photos in Figure XII 
indicates a typical turbulent flow pattern. As speed in- 
creased the mere fully turbulent flow region was approached 
with the added increased effect of the trip wire with speed 
thus causing the point of separation to move aft on the body, 
as it does in a typical transition from laminar to turbulent 
flow. The two theoretical separation criteria are insensi- 
tive to changes in speed in the range investigated as has 
been previously discussed. 

The recommendation here is obvious; the experiments 
must be conducted at higher speeds to further insure turbu- 
lent flow. In the present experiments the speeds obtainable 
were directly related to the ability to observe the dye flow, 
the dye being dissipated with increased speed and almost 
wholly indistinguishable beyond 1.75 knots. The speed limi- 
tation was due in part also to large carriage vibrations at 
higher speeds than those used. With a more refined method of 
turbulence stimulation and dye=>flow observation techniques 
higher speeds could be attained leading to perhaps more con- 
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elusive and plausible results relating the effect of changing 
speed on the separation point than those herein achieved. 

Corrected Body Profiles and Pressure-Velocity Profiles 

Equivalent or corrected body profiles shown in Figures 
XVII through XXII are the results of Computer Program II. 

From the sample of these results given in Table G-I in 
Appendix G it can be noted that there is only a very slight 
variation in corrected profile on a given body with varia- 
tion velocity. The corrected body profile radii decreased 
slightly with increase in speed. The variation is considered 
to be beyond the accuracy of the model construction and the 
full scale drawings so the corrected body profile was plotted 
the same for all velocities on a given body. In Figures XVIII 
through XXII the corrected body is plotted beyond the separa- 
tion point however, this in fact does not exist since this 
profile is based upon modified potential flow theory which 
does not hold after separation occurs. 

The pressure-velocity profiles for the bodies shown 
in Figures XIII through XVII are the results of Computer Pro- 
gram III. From the sample of these results given in Table G-I I 

U/ 

in Appendix G only a slight variation in Cp and Uoc occurs 
with variation of Uoo or model velocity. This is a result of 
the aforementioned slight variation in body profile with model 
velocity since the output of Program II, Y - NORM, is the in- 
put to Program III. Again this slight variation is considered 
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beyond the accuracy of the model drawings and certainly is 
not detectable when plotted at the scales used in Figure XIII 
through XVII. Thus the profiles appear the same for all ve- 
locities on a given body. 

Since neither depth of submergence nor hydrostatic 
pressure enter into the calculations of the corrected body 
profile or pressure- velocity profiles none of the figures 
show a variation with depth. 

In summary, the experimental results indicate that 
separation appears to begin at a critical angle of about 
25 degrees irrespective of body streamlining. The correla- 
tion of separation with speed and depth is not clearly seen 
since the speed and depth range should be larger for more 
conclusive results. The results as pointed out in the in- 
troduction are qualitative in nature and are not be be con- 
strued as quantitative evidence. Such evidence must be pro- 
vided in later works. The primary recommendation is that 
further studies be carried out at greater speeds with im- 
provements in the technique of observing the separation 
point. Stratford’s separation criterion, which is more ri- 
gorous than Buri’s appears to give good agreement with the 
experimental observations, but it should be further evaluated 
by using experimentally determined pressure-velocity dis- 
tributions. 
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V - CONCLUSIONS 



The conclusions which were arrived at in IV are summa- 
rized as follows: 

1. The observed position of separation was determined 
within an accuracy of - 1/2 inch on the body and the body 
angle at this position was thereby determined to - 3? 

2. Initial separation occurs at body angles of be- 
tween 23° and 28° regardless of the amount of body streamlining. 

3. Buri’s criterion for separation consistently predicts 
separation well forward of the observed separation point and 
thus does not compare favorably with the observed results. 

4. Stratford’s criterion evaluated from the calcu- 
lated pressure-velocity distributions over the bodies gives 
reasonably close agreement with the location of the observed 
separation point. 

5. Results of analyzing the separation position lo- 
cation show that it moves aft with increase in speed but this 
may be due to the fact that at the speeds investigated the 
Reynolds numbers were in the transition region between laminar 
and turbulent flow. 

6. The theoretical criteria evaluated are both in- 
sensitive to body speed changes and are independent of the 
depth of submergence of the body. 
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7. The observed separation point does not move in a 
consistent manner on the body with a variation in depth of 
submergence of the body. 

8. The corrections applied to the body surface to 
account for boundary layer effects showed little change with 
change in body speed for a given body. 

9. The pressure-velocity profiles over the bodies 
calculated from modified potential flow showed little change 
with change in body speed for a given body. 
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VI - RECOMMENDATIONS 



The recommendations made in IV are summarized as 
follows: 

1. Further experimentation should employ a larger 
model if possible and an improved means of observing the 
flow about the separation point. 

2. A greater range of speeds and depths should be 
investigated in future work to increase the reliability of 
the data obtained. Higher speeds consistent with the a- 
bility to observe the dye-flow are necessary to insure that 
flow is definitely in the turbulent region. 

3. If higher speeds can not be attained improve- 
ment should be made in the method of inducing turbulence 
through stimulation. 

4. If dye studies are continued a minimum of 4ye 
should be used to prevent dye outside the boundary layer 
from obscuring the observations. 

5. Pressure measurements should be made in and 
ahead of the separation zone to experimentally determine 
the pressure-velocity distribution over the bodies and also 
serve as another method of locating the separation point. 
This could be compared with the results of this work and 
used to re-evaluate the separation criteria on the basis of 
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the actual pressure-velocity profiles rather than on calcu- 
lated profiles. This appears to be the next logical step 
as the initial determination of the separation zone has been 
made in this work. 

6. Recommendations relating to drag force measure- 

i 

ments are included in Appendix E. 
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APPENDIX A 

OFFSETS AND OTHER CHARACTERISTICS FOR MODEL AND STERN SECTIONS 



TABLE A -I - OFFSETS 



Body: 


Nose Section 


Stern 1 


Stern 2 


Stern 3 


Stern 4 


Stern 5 


St at ion 


Spacing: 0.50" 


0.50" 


0.50" 


0.50" 


0.50" 


0.50" 


Station 


0 at: Nose 


Tail 


Tail 


Tail 


Tail 


Tail 


Station: 


Radius ( Inches ) 


Radius 


Radius 


Radius 


Radius 


Radius 




0 


0.00 




0.00 


0.00 


0.00 


0.00 


0.5 


0.43 




0.80 


0.65 


0.63 


0.32 


1 


0.68 


CO 

Uj 

X 


1.05 


0.87 


0.85 


0.47 


1.5 


0.83 


1.23 


1.03 


0.98 


0.57 


2 


0.94 


o 

> 


1.37 


1.16 


1.09 


0.69 


3 


1.12 




1.58 


1.36 


1.27 


0.87 


4 


1.27 


o 


1.72 


1.52 


1.41 


1.00 


5 


1.37 


Q 

CVJ 

II 


1.84 


1.66 


1.53 


1.12 


6 


1.47 


1.93 


1.77 


1.63 


1.23 


7 


1.54 


1.98 


1.84 


1.72 


1.33 


7.26 


— 


lo 

ID 


2.00 


— 


— 


— 


8 


1.62 


o 




1.91 


1.78 


1.42 


9 


1.67 


< 




1.95 


1.84 


1.50 


10 


1.72 


oc 




1.98 


1.89 


1.57 


11 


1.76 


li. 




1.99 


1.93 


1.63 


11.30 


— 


0 




2.00 


— 


— 


12 


1.80 








1.96 


1.69 


13 


1.83 


UJ 

CY 






1.98 


1.74 


14 


1.87 


Ui 






1.99 


1.79 


14.78 


— — 


X 

n 






2.00 


— 


15 


1.90 


LL 

l/) 








1.83 


16 


1.93 


r 








1.87 


17 


1.95 


lu 








1.90 


18 


1.96 


X 








1.92 


19 


1.97 










1.94 


20 


1.99 










1.96 


21 


2.00 










1.98 


22 


— 










1.99 


22.78 


— 










2.00 
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TABLE A-II- OTHER CHARACTERISTICS 



Nose/Tail 

Radius 

(inches) 


Length 

(inches) 


Wetted Surface 
Depths 1 and 
(sq.ft . ) 


Wetted Surface 
2 Depth 3 

(sq.ft. ) 


Nose 0.43 


10.5 


0.726 


0.640 


Cylindrical Section — 


19.0 


1.653 


1.311 


Stern 1 2.00 


2.00 


0.175 


0,158 


Stern 2 1.47 


3.63 


0.245 


0.219 


Stern 3 0.98 


5.65 


0.387 


0.340 


Stern 4 0.99 


7.39 


0.513 


0.454 


Stern 5 0.33 


11.39 


0.741 


0.664 


Body 1 


31.50 


2.554 


2.-109 


Body 2 — 


33.13 


2.624 


2,170 


Body 3 


35.15 


2.766 


2.291 


Body 4 — 


36.89 


2.892 


2.405 


Body 5 


40.89 


3.120 


2.615 


Body means total of nose, 


cylindrical section, and 


stern shape. 


Wetted surfaces submerged 


(depths 1 


and 2) were calculated 


using ten-ordinate Simpson's rule with half-stations on 


body perimeters 


Spacings of 1 inch were used on nose 


> and stern 5, 


3/4 inch on 



stern 4, 1/2 inch on sterns 2 and 3. Stern 1 was calculated 



directly since it was a hemisphere. 

Surfaced wetted surfaces (Depth 3) were calculated by computing 
non-wetted surface Using trapezoidal rule on non-wetted arc lengths 
and subtracting from submerged wetted surface values. 

These calculations are not shown in this report. 
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TABLE A-III 



GAUSSIAN STATION, RADIUS, AND BODY ANGLE 



Gaussian X Body i Body 2 Body 3 



Station L 


X 


Radius 


: Angle X 


Radius 


Angle 


X Radius 


Angli 


-.990 


.005 


0.17 


0.37 


49.5 


0.18 


0.37 


49.4 


0.19 


0.40 


47.3 


-.944 


.028 


0.87 


0.88 


26.0 


0.92 


0.91 


25.0 


0.97 


0.94 


23.9 


-.866 


.067 


2.12 


1.30 


14,2 


2.23 


1.31 


12.7 


2.36 


1.37 


11.7 


-.756 


.122 


3.85 


1.60 


7.5 


4.05 


1.63 


7.1 


4.30 


1.67 


6.3 


-.618 


.191 


6.02 


1.81 


4.5 


6.33 


1.84 


4.4 


6.72 


1.87 


3.6 


-.458 


.271 


8.54 


1.96 


2.6 


8.98 


1.98 


1.8 


9.52 


1.97 


1.5 


-.282 


.359 


11.31 


2.00 


0.0 


11.90 


2.00 


0.0 


12,62 


2.00 


0.0 


-.096 


.452 


14.25 


2.00 


0.0 


14.99 


2.00 


0.0 


15.90 


2.00 


0.0 


+ .096 


.548 


17.25 


2.00 


0.0 


18.14 


2.00 


0.0 


19.24 


2.00 


0.0 


+ .282 


.641 


20.19 


2.00 


0.0 


21.23 


2.00 


0.0 


22.52 


2.00 


0.0 


+ .458 


.729 


22.96 


2.00 


0.0 


24.15 


2.00 


0.0 


25.62 


2.00 


0.0 


+ .618 


.809 


25.48 


2.00 


0.0 


26.80 


2.00 


0.0 


28.43 


2.00 


0.0 


+ .756 


.878 


27.65 


2.00 


0.0 


29.08 


2.00 


0.0 


30.85 


1.94 


5.1 


+ .866 


.933 


29.38 


2.00 


0.0 


30.90 


1.78 


13.1 


32.7$ 


1.62 


15.5 


+ .944 


.972 


30.63 


1.63 


35.0 


32.21 


1.33 


27.8 


34.18 


1.15 


24.5 


+ .990 


.995 


31.33 


0.78 


67.2 


32.95 


0.69 


59.6 


34.96 


0.57 


51.5 



Gaussian X Body 4 Body 5 



Station 


L 


X Radius 


Angle 


X Radius 


Angle 


-.990 


.005 


0.20 0.41 


47.1 


0.21 0.42 


46.5 


-.944 


.028 


1.02 0.97 


22.7 


1.13 1.00 


21.6 . „ 
10.0 * orB 


-.866 


.067 


2.48 1.38 


11.5 


2.75 1.42 



-.756 


.122 


4.51 


1.67 


5.7 


5.00 


1.72 


5.1 


X is inches from 


-.618 


,191 


7.05 


1.88 


3.4 


7,81 


1.91 


2.7 


nose of body. 


-.458 


.271 


10.00 


1.99 


1.0 


11.08 


2.00 


0.0 


Radius is in 


-.282 


.359 


18y25 


2.00 


0.0 


14.69 


2.00 


0.0 


inches. 


-.696 


.452 


16.6$ 


2.00 


0.0 


18.50 


2.00 


0.0 


Angle is in 


+ .096 


.548 


20.20 


2.00 


0.0 


22.39 


2.00 


0.0 


degrees. 


+ .282 


• .641 


23.64 


2.00 


0.0 


26.20 


2.00 


0.0 




+.458 


.729 


26.89 


2.00 


0.0 


29.81 


1.98 


1.0 




+.618 


.809 


29.84 


1.99 


1.0 


33.07 


1.85 


4.1 




+ .756 


.878 


32.38 


1.84 


6.2 


35.89 


1.57 


8.3 




+ .866 


.933 


34.41 


l.£2 


12.0 


38.14 


1.18 


'12.9 




+ .944 


.972 


35.87 


1.10 


22.0 


39,76 


0.73 


20.4 




+ .990 


.995 


36.69 


0.55 


53.0 


40.67 


0.28 


36.4 
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APPENDIX B - DESCRIPTION OF TOWING TANK AND INSTRUMENTATION 



An early description of the Ship Model Towing Tank, 
prior to the installation of the towing carriage, is given in 
a paper entitled "The Ship Model Towing Tank at MJ.T.," Society 
of Naval Architects and Marine Engineers, 1953. This paper 
does, however, give an up to date description of the tank di- 
mensions and the speed measuring and counting equipment. Since 
the towing carriage was installed late in 1961, a brief de- 
scription will be included here. 

The carriage is constructed of aluminum and rides on 
two hollow cylindrical stainless steel rails, one overhead and 
one mounted on the tank’s side wall. The overhead rail is 
simply for support and guidance of the carriage which rides it 
on six fibreglass wheels, two above and one below the rail at 
either end of the carriage frame. The carriage is also sup- 
ported by the side rail with three wheels riding on the after 
end and a soft rubber drive wheel connected via a coupling to 
an electric motor on the forward end. The electric drive motor 
is controlled by the console at one end of the tank via cables 
attached to pulleys which slide along a smaller rail mounted 
on the tank wall above the side rail. The control cables are 
free riding and consequently as the carriage moves from one 
end of the tank to the other the speed decreases as more cable 
length is either drawn, out from one end or piled up. A dif- 
ferential pulley system was installed which helped in alleviating 
this situation to some degree. 
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The carriage speed can be controlled either manually by 
setting a dial on the console or automatically at a set speed 
or to position the carriage over a self propelled model. The 
automatic control was inoperative during the period of the ex- 
periments included in this thesis, and speed was controlled 
manually. The speed measuring apparatus includes a type 
500 Hewlett Packard frequency meter and a type 521 Hewlett 
Packard counter described in greater detail in the foremen- 
tioned reference- 

Figure B-I shows a schematic drawing of the model dye 
system which is self-explanatory and figure B-II shows the 
details of the drag measuring towing dynamometer. 
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FIGURE B-T , MODEL DYE SYSTEM 
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FIGURE B-2E DYNAMOMETER DETAILS 
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APPENDIX C - COMPUTER PROGRAMS 



The IBM 7090 electronic computer located at the M.I.T. 
Computation Center was utilized as an aid in reducing drag 
data and in computing the theoretical portions of the thesis. 

Three programs were written: 

Program I - Reduction of drag data. 

Program II- Calculation of turbulent boundary layer charac- 
teristics on the bodies tested. 

Program III- Calculation of the pressure velocity distribu- 
tion on the bodies tested. 

These programs were written in 7090 FORTRAN. On the 
next several pages the following are given for each of the 
three programs: 

1. FORTRAN program listing. 

2. A list of program symbols and their meanings. 

3. A logical sequence flow chart of the program. 

Both FORTRAN and compiled binary machine language card 
decks and listings are available for these programs in the 
Department of Naval Architecture and Marine Engineering 
computer program files. 
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PROGRAM I - LISTING OF SYMBOLS 



SYMBOL 

RO 

VIS 

SLEN(N) 
WSUR(N ,M) 

DT(L,M,N) 

V' (L,M,N) 
N 

M 

L 

SL 

DSTR 

DTC 

SLR 

RM 

A 

CT 

CFSG 

CFS 

CR 



MEANING 



UNITS 



Mass density of tank water 
Kinematic viscosity of tank water 
Body length 

Wetted surface of body at 
particular depth 

Measured drag 

Model velocity 

Body index 

Depth index 
Speed index 

Lenth of strut submerged 
Strut drag 
Corrected body drag 
Speed - length ratio (V/Y~L) 



slug/ft. 3 

ft. 3 /sec. 

ft. 

ft. 2 

lbs. 

kts. 



inches 



lb. 



lb. 






Reynolds number on SLEN (N) 

1/2 $ S V 2 lb. 



Total drag coefficient 

Friction drag coefficient using 
Shultz-Grunow . formulation 

Friction drag coefficient using 
Schoenherr f ormul at ion 

Residual drag coefficient 
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LOGICAL SEQUENCE FOR PROGRAM 1 



START 



I 



READ f , (RO) 
and if , (VIS) 

f I 



READ LENGTH OF 
BODY, (SLEN(N)) 



READ VELOCITY 
V(L,M,N) 



SET 


BODY INDEX 




N-l 



d> 



PRINT BODY 
NUMBER - N 



I 



SET DEPTH INDEX 
M-l 






PRINT DEPTH- M 
and OUTPUT HEADING 




IF M-2 




© 

i 



COMPUTE STRUT DRAG, 
(DSTR) , Equation (27) 



DSTR-0.0 



COMPUTE CORR. 
DRAG, (DTC) 



1 


l 




Y. . 

COMPUTE V/L^CSLR) 
and REYNOLDS NUMBER, (RN) 


READ MEASURED 

TMJ bn. TYT ( T M 











COMPUTE C t ,(CT) 



COMPUTE Cf from Schultz- 
Grunow Formula, (CFSG) 

Equation (30) 



CFSG-CFS 



COMPUTE Cf from Schoenherr 
Formula using CFSG 
Equation (Z“5) 



^^J ^7(ABS.VAL. (CFSG-CFS) -.OOl ) 




STRUT LENGTH, 


1 ' 


STRUT LENGTH, 




fl>RINT V(L,M,N) ,SLR,DT(L,M,N) ,DTC, 


SL— 0.0 


. 


SL— 6 . 0 




CT ,CF,CR, and RN 



V 



STRUT LENGTH, 
SL-14.0 



( IS I^6> ><N^>» f^°Ln 



SET VELOCITY INDEX 
L=1 





(!SM-3>X@Hm^ 




I 

© 



N-N+lH — <^x js~ N-"5~> »<jE^>END 
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PROGRAM I FOR CALCULATION OF CR ON BODY 



FORTRAN 

DIMENSION DT(6,3,5),V(6,3,5),SLEN(5),WSUR(5,3) 

READ 24 , RO , VIS 

24 FORMAT (F5.3, E10.4) 

READ 25, (SLEN(N),N-1,5) 

25 FORMAT (5F4.2) 

READ 26, ( (WSUR(N,M) , N-l ,5) , M-l ,3) 

26 FORMAT (14F5.3) 

READ 27, ( ( (DT(L,M,N) , L-1,6) M-l ,3) , N-l , 5) 

27 FORMAT (I8F4.2) 

READ 28, ( ( (V(L,M,N) , L-1,6), M=1 , 3) ,N=1 , 5) 

28 FORMAT (18F4.2) 

DD 50 N-l, 5 
PRINT 32, N 

32 FORMAT (14H1 BODY NUMBER-11) 

DP 50 M-l, 3 

PRINT 33, M 

33 FORMAT (8H0 DEPTH-11) 

PRINT 34 

34 FORMAT (111H0 VEL. (KT) REYNOLDS NO. SPEED/LENGTH TOT.DRA 

1 6 TOT. DRAG (CORR) CT CF CR/) 

IF (M-2) 10,11,12 

10 SL-14.0 
GO TO 14 

11 SL-6.0 
GO TO 14 

12 SL-0.0 

14 DO50 L-1,6 

IF (M-2) 61,61,60 

60 DSTR-0.0 
GO To 62 

61 DSTR- (RO* ( (1 .689*T(L,M,N) )**2. 0)* ( (0. 064*SL)+1 . 17) )/288 . 0 

62 DTC-DT (L , M , N) - DSTR 
SLR-V(L,M,N)/SQRTF(SLEN (N) ) 

RN— 1 . 689* V (L^M,N)*SLEN (N)/VIS 

A-0.5*RO*WSUR(N,M)*(1.689*V(L,M,N))**2.0 

CT-DTC/A 

Xl-0 . 434294*LOGF (RN ) 

CFSG— 0. 0725/ (Xl-2. 0)**2. 0 

13 X2— 0. 434294*LOGF(RN*CFSG) 

CFS-0.0585/x2**2. 0 

IF ( ABSF (CFSG- (0 . 001E-03 ) ) 7 , 7 , 1 5 

15 CFSG-CFS 
GO TO 13 

7 CR-CT-CFS 

PRINT 35, V(L,M,N) , RN,SLR,DT(L,M,N) ,DTC,CT ,CFS ,CR 

35 FORMAT (3H F4.2,7H 1PE10.3,9H 0PF5.2,1CH 

1 0PF5 . 3 , 8H 0PF3.3 , 10H 1PE10.3,4H 1PE10.3,5H 

21PE10.3) 

50 CONTINUE 
CALL EXIT 
END 
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PROGRAM II - LISTING OF SYMBOLS 



Symbol 


Meaning 




Units 


P 


Reciprocal of power in 
velocity profile 




- 


VIS 


Kinematic viscosity of tank water 


2 

ft. /sec. 


SLEN(N) 


Body length 




ft* 


V(L,M,N) 


Model velocity 




kts* 


X(I,N) 


Distance from body nose 




inches 


R(I,N) 


Radius of body at X(I,N) 




inches 


AL ( I , N ) 


Cosine of angle between body sur- 
face and body axis at X(I,N) 


- 


N 


Body index 




- 


M 


Depth index 




- 


L 


Speed index 




- 


I 


Location on body index 




- 


RN 


Reynolds number on X 




- 


PBLT 


Boundary layer thickness 
flat plate ( opl) 


on 


inches 


PDT 


Displacement thickness on 
plate (£* pi) 


flat 


inches 


BLTO 


Estimate of boundary layer 
thickness on body 


inches 


BLT 


Boundary layer thickness 
body ( o ) 


on 


inches 


A 


Ratio of boundary layer thickness 
to model radius 


— 


DLNP 


An pi 




— 


DLTNP 


AZn p i 






DT 


Displacement thickness on 
(£*) 


body 


inches 
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PROGRAM II (CONT * D) 



Symbol 


Meaning 


Units 


OMNP 


n 

Pi 


— 


OMTNP 


2n p i 


— 


BMT 


Momentum thiickness on body (0) 


inches 


RC 


Body radius corrected to 
account for displacement 
thickness 


inches 


Y 


Normalized radius 
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LOGICAL SEQUENCE FOR PROGRAM II 



START 

_X_ 

READ 1/n , (P) ; if, (VIS ) ; and 
BODY LENGTH, (SLEN(N)) 



© 



3 



COMPUTE © ON FLAT PLATE , 
(BMT) , Equation (3a) 



E 



READ VELOCITY 




SET £ ON BODY* £ ON FLAT PLATE 


V(L,M,N) 




(BLTO-PBLT) 



IE AD LOCATION ON BODY 
X(I,N) 

i ; — 

READ BODY RADIUS 
R(I,N) 



COMPUTE £ ON BODY 
(BLT), Equation (l) 



BLTO-BLT 



<^><(^IF~ (ABS.VAL. (BLTO-BLT)-.OOl/) 



READ COSINE OF ANGLE BETWEEN] 
BODY SURFACE AND BODY AXIS 
AL (N , I ) 




or 



SET BODY INDEX 

_N=1. 



N-N+l h- 

SET DEPTH INDEX 
M-l 

* — | m-m+i |< - 



<D 



COMPUTE £ * ON BODY , 
(DT) , Equation (2) 



COMPUTE © ON BODY, 
(BMT), Equation (3) 



SET VELOCITY INDEX 
L-l 



L-Lfl 



~® 


COMPUTE CORRECTED BODY RADIUS, 
(RC), Equation (4) 




'C 



IMPRINT OUTPUT HEADINGS 



SET LOCATION INDEX 
1- 1 



I-I+l 



© 



COMPUTE REYNOLD’S 
NUMBER, (RN) 



COMPUTE NORMALIZED BODY RADIUS, 
(Y), Equation (4a) 






T 



PRINT X(I,N),R(I,N)„RC, BLT , 
DT ,BMT , RN , and Y 



(IS L-3> <jfE^>- (lS 1-16) 



COMPUTE £ ON FLAT PLATE > 
(PBLT), Equation (Id ) 



I 




COMPUTE <£ * ON FLAT PLATE 
(PDT) , Equation (2b) 



Gs M-3 ) 

W <^> 

© <» N-5 ~> K@>— 




es; 



END 
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PROGRAM II, CALCULATION OF BODY DISPLACEMENT THICKNESS , ETC. 



FORTRAN 

DIMENSION V(3 , 3 , 5) ,X(16,5) ,R(l6 , 5) ,SLEN(5) ,AL(l6,5) 

READ 21 ,P,VIS , (SLEN (N) ,N-1 , 5) 

21 FORMAT (F3 . 1 ,E10. 4 , 5F4. 2) 

READ 22, (((V(L,M,N),L“1,3),M-*1,3),N-1,5) 

22 FORMAT (18F4.2) 

READ 23, ( (X (I, N), 1-1,16) ,N-1 ,5) 

23 FORMAT (14F5. 2) 

READ 24 , ((R(I,M),I— 1,16), N- 1,5) 

24 FORMAT (16F4.2) 

READ 25, ((AL(I,N),I-1,16),N-1,5) 

25 FORMAT (14F5.3) 

DO 50 11-1,5 

DO 50 M-1,3 
DO 50 L-1,3 
PRINT 100, N 

100 FORMAT (14H2 BODY NUMBER-11) 

PRINT 110, M 

110 FORMAT (13H0 BODY DEPTH-11) 

PRINT 101 , V (L,M,N) 

101 FORMAT (16H0 VELOCITY (KTS )-F4 , 2 ') 

PRINT 102 

102 FORMAT (95H0 X(INCH) RADIUS CORR. RADIUS DELTA DELTA 
1STAR THETA REYNOLDS NO, Y-NORM/) 

DO 50 1-1,16 

RN-1.689*V(L,M,N)*X(I,N)/(VIS*12.0) 

PBLT-0.38*X(I,N)/RN**0.2 
PDT-PBLT/ (l , 0+P) 

BLTO-PBLT 

12 BLT=PBLT/(l.0+(BLTO/(3.0*R(I,N))))**0.8 
IF (ABSF(BLTO-BLT)-. 001) 7,7,10 
10 BLTO-BLT 
GO TO 12 
7 A-BLT/R ( I , N ) 

DLNP-PDT/PBLT 

DLTNP-DLNP* (P+1 , 0)/ (2. 0 s * P+1 . 0) 

DT-PDT* (1 , 0+A-SQRTF(l , 0+2. 0* A* (l . 0-DLNP)+ (A* *2 . 0)* (1 . 0-DLTNP) ) )/ (A 
1 *DLNP* (l . 0+A/3 . 0)**0„ 8) 

OMNP-p/ ( (P+1 .0)* (P+2 . 0) ) 

OMTNP-OMNP* (P+2 . 0 ) /2 . 0*P+ 1.0) 

BMT-BLT* ( OMNP+ A* OMTNP/ 2 . 0 ) 

RC«R(I,N)+AL(I,N)*DT 
Y-RC/ (SLEN (N) *6 . 0 ) 

PRINT 103,X(I,N),R(I,N) , RC , BLT , DT , BMT , RN , Y 

103 FORMAT (3H F5.2,6H F4.2,7H F5.3,9H F5.4,7H 

1 F5.4 ,6H F5.4,5H 1PE10.3,4H 0PE12.5) 

50 CONTINUE 
CALL EXIT 
END 



71 





PROGRAM III - LIST OF SYMBOLS 


Symbol 


Meaning Units 


(I)+(J) 

X(I) 


Gaussian abscissae, number (1-16) 

Normalized Gaussian abscissae based on model 
length of two, from -1 to+ 1 


Y(I) 


Normalized ordinate of body meridian profile 


R(I) 


Gaussian weighting functions 


RO 


Water density slug/ ft** 


PATD(M) 


Static pressure (gage) PsiG 


V(L,M,N) 


Absolute model velocity kts 


N 


Body number 


M 


Depth 


L 


Speed 


AL(I) 


Cosine of angle between tangent to body profile 
and x-axis 


A(I) 


Ml) ratio of length to diameter 


B(I) 


k 1 (I), longitudinal virtual mass coefficient 


G(J,I) 


g(x,t), equation of prolate spheroid whose ends 
coincide with ends of the given body, and which 
intersects body at x -t. 


C(J,I) 


K (x,t), matrix of functions of body 


D(J,tt) 


K / (x,t) , Matrix of functions of prolate spheroid 


SUM(I) 


£( K-K») 


E(K, I) 


E(n,t), error functions 


H(I) 


k t + E 
1 + kj (x) 


K 


Iteration number 
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4 

I I 




■ 




PROGRAM III - (CONT’D) 



Symbol 


Meaning 


Units 


SUM (I) 


£(En . - E ni ) 
nj ni 




0(1) 


^Eni 




U(D 


, ratio of flow velocity over 
stream velocity. 


body to free 


p(i) 


p/q j pressure distribution 




VEL(I) 


Absolute flow velocity over body 


kts . 


PCC(I) 


Gage pressure on body 


lb/ft 2 
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LOGICAL SEQUENCE FOR PROGRAM III 



Start 





Read X(I), R(I), RO, 
PATD(M), V(L,M»N) 


0 


* 5 , 


r 




Read N,M,L, AL(17) ' 




3 


f 




Read Y(I),AL(I) ( 



Q 




Is J - 16 



3 



r-0 -4^ 



jHhl 




Compu t e ^ (K-K * ) 



P - 4 



0 H 




Compute E , 

E(l,T) eq. (19) 






P-P+l 



I -1 



0 - 



Compute H(I) - 
B(17) + E(l,I) 



i + fi. Ci73 



-0 




Is I - 16 





na 



I-I+l 




9 



K - 1 





Compute ^K(E nS - E n ^ 



(lK(J,I)x(E(K,J)-E(K ,1))) 



Compute X(x)^(A(I)) 
f 4 W equation (17) 



Compute K^(x),(B(I)) 
( 5 equation (18) 

V _ ' Compute B(17) 



-0 



-( 7 ) 



( 



Is J - 16 




) 




J - J-ll Is I - 16 ) 



9B 



<0> 



(J) 



Compute g (x , t ) , 
(G(J,I)) eq. (16) 



Compute K(x,t), 
(C(J 9 1) eq , (14) 



\ 


f 


Compute E n . t j 
(E(K+ l f I))“ 1,1 
eq. (20) 


' 


f 






I 1 - 


I 41 











fis k -5 3 





! p 


> 


t 


K- 


K+l 


I- 


1 




Comput e K * (x , t ) , 
(D(J , I) ) eq. (15) 




13 



9A 




(eoEtinued on next page) 
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Print N,M, L and 
output heading 



0 



Compute E i 
E(K,I),K« 1..5 



Compute U(I) 

-d+ki + ZE ni ) COS 
-(1+B(17)+2E(K,I))AL(I) 



Compute p/q 
P(I)-1-U(I)2 



Compute v bs 
VEL(l)-U(l)xV{L s M f N) 



Compute P„,_,. 

* ^d A<5 E q 

PCC(I)-P(I)xROxV(L s M»N)^ 

2 

+ PATD(M) 



Print X(I),Y(I),U(I) 
P(I),PCC(I),VEL(I) 
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PROGRAM III CALCULATION OF PRESSURE DISTRIBUTION 



FORTRAN 

DIMENSION X(16),R( 16) , Y(16) , AL(16 , 5) , A(17) .Dim) ,C(16 , 16) , D(16 , 16) 
1,SUM(16),E(5,16),H(16),0(16),U(16),S(16) ,P (16) ,G(16 , 16) ,PATD(3 ) , 
2V(3 ,3 , 5) , VEL(16) ,PCC (16) 

READ l,(X(l),R(l),T-l,16) 

1 FORMAT (6E12.5) 

READ2, RO , (PATD (M) , M-l , 3 ) 

2 FORMAT (F5. 3, 3E1 1.4) 

READ3, (((V(L,M,N),L-13),M-1,3),N-1,5) 

3 FORMAT (18F4.2) 

READ7, ((AL(I,N),I-1,16),N-1,5) 

7 FORMAT (14F5. 3) 

4 READ5,N 9 M s L 9 A(l7) 

5 FORMAT (311, El 2. 5) 

READ6, (Y(I), 1-1,16) 

6 FORMAT (6E12. 5) 

DO 10 1-1,16 

A(D-SQRTFC(X(I)+I.)*(l.-X<r))/(Y(I)**2.)) 

10 CONTINUE 
DO 20 1-1,17 

B(D-(A(r)*LOGF(A(T)+SQRTF(A(l)** 2.-1 . ) )-SQRTF(A (l)**2.-l .))/ 

1 ( (A(l)**2. )*SQRTF(A (l)**2. -1 . )-A(r)*LOGF(A'(l)+SQRTF(A (l)**2. -1 . ) ) ) 
20 CONTINUE 

DO 30 1-1,16 
DO 30 J-1,16 

G(J„I)-(Y(T)**2.)*(X(J)+1.)*(1.-X(J))/((X(T)+1.)*(1.-X(I))) 

C(J „T)-R(J)* (Y(J)**2. )/( ( (X(J)-X(l) )**2.+Y (J)**2. )**1 . 5) 

D(J j,T)-R(J)*G(J ,T)/( ( (X(J)-X(I) )**2.+G(J 9 D )**1 . 5) 

30 CONTINUE 

DO 40 1-1,16 
SUM (I)- 0.0 
DO .39 J-1,16 
S(J)-C(J,I)-D(J,I) 

SUM(r)-SUM(l)+S (J) 

39 CONTINUE 
K-l 

E(K,r)-l.-(((L+B(l7))/2.)*SUM(I)Hl.4B(l7))/(l.+B(I)) 

H(l)-(E(K 9 r)+B(17))/(l.+B(17)) 

40 CONTINUE 

DO 50 K-l, 4 
DO 50 1-1,16 
SUM(1)=0. 0 
DO 49 J-1,16 

S (J)— C (J jT)* (E (K, J)-E(K,I) ) 

SUM(I)- SUM (I) + S(J) 

49 CONTINUE 
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PROGRAM III (CONT"D) 



E <Kfl,I)-H (X) *E (K ,1) - . 5*SUM (r) 

50 CONTINUE 
PRINT 52, N 

52 FORMAT (14H2 BODY NUMBER-11) 

PRINT53 U 

53 FORMAT (8H0 DEPTH-Il) 

PRINT54 , V (L,M,N) 

54 FORMAT (8H0 SPEED- F4. 2) 

PRINT55 

55 FORMAT (83H0 X(I) Y(I) U(I) P (X) 

1 PCC(I) VEL(X)/) 

DO 60 1-1,16 

0(l)=E(l,I)+E(2,I)+E(3,I)+E(4,I)+E(5,I) 

U(l)-(l.+B(l7)+0(l))*AL(I,N) 

P(l)-1.-(U(I))**2. 

VEL(l)— U(I)*V (L,M,N) * 

PCC (T)-P (!)♦(. 5*RO* ( (V (L,M,N)*1 .689)**2. ) )+PATD(M)*144. 

57 PRINT 59, (X (I ) , Y (I ) , U (I ) , P (X) , PCC (I ) , VEL (T ) ) 

59 F0RMAT(6E14. 5) 

60 CONTINUE 

GO TO (4, 4, 4, 4, 61), N 

61 GO TO (4,4,62) ,M 

62 GO TO (4,4,63) ,L 

63 CALL EXIT 
END 
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APPENDIX D - EXPERIMENTAL DATA 
TABLE D-I _ EXPERIMENTAL DRAG DATA 



Dial Speed Voltmeter Voltmeter Voltmeter Measured 
Run Stern Depth Speed (KTS) Zero(V) AVG (V) Dlif . (V) Drag(lbs) 



1 


3 


1 


20 


0.80 


- 43.7 


- 26.4 


17.3 


0.21 


1 A 


/ 








20 


0.77 


- 45.0 


- 27.9 


17.1 


0.21 


2 










24 


1.12 


- 46.5 


- 20.4 


26.1 


0.33 


2 A 










24 


1.12 


- 45.7 


- 18.5 


27,2 


0.34 


3 










28 


1.48 


- 45.0 


- 14.75 


30.25 


0.38 


3 A 










28 


1.48 


- 43.5 


- 11.7 


31.8 


0.40 


4 






2 


20 


0.75 


- 45.1 


- 29.2 


15.9 


0.20 


4 A 


3 






20 


0.74 


- 47.7 


- 30.8 


16.9 


0.21 


5 










24 


1.09 


- 45/5 


- 22.6 


22.9 


0.24 


5 A 










24 


1.09 


- 44.7 


- 18.5 


26.2 


0.33 


6 








< 


28 


1.50 


- 43.0 


- 18.4 


24.6 


0.31 


6 A 






i 




28 


1.51 


- 43.5 


- 19.0 


24.5 


0.31 


7 






r 


; 


20 


0.77 


- 46.0 


- 30.7 


15.3 


0.19 


7 A 










20 


0.79 


- 46.1 


- 30.9 


15.2 


0.19 


8 










24 


1.16 


- 45.5 


- 23.6 


21.9 


0.27 


8 A 










24 


1.17 


- 45.0 


- 25.1 


19.9 


0.25 


9 


1 


t 






28 


1.62 


- 46.0 


- 19.6 


26.4 


0.33 


9 A 


3 


r 


i 


28 


1.58 


- 45.0 


- 21.0 


24.0 


0.30 


10 






r 




20 


0.77 


- 45.0 


- 23.2 


21.8 


0.28 


10 A 










20 


0.77 


- 46.0 


- 21.2 


25.0 


0.31 


11 A 










24 


1.16 


- 46.0 


- 24.0 


22.0 


0.28 


11 B 










24 


1.21 


- 46.5 


- 25.8 


20.7 


0.26 


12 










28 


1.58 


- 46.0 


- 19.7 


26.3 


0.33 


12 A 






r 


i 


28 


1.58 


- 45.5 


- 20.0 


25.5 


0.32 


13 






2 


20 


0.75 


- 35.0 


- 22.6 


12.4 


0.16 


13 A 










20 


0.75 


- 36.0 


- 24.2 


11.8 


0.15 


14 


1 






24 


1.19 


- 35.0 


- 13.2 


21,8 


0.28 


14 A 










24 


1.11 


- 37.5 


- 15.2 


22.3 


0.28 


15 






\ 


t 


28 


1.61 


- 42.5 


- 12,5 


30.0 


0.38 


15 A 






2 


28 


1.61 


- 44.0 


- 12.6 


31.4 


0.39 


16 






] 


L 


20 


0.74 


- 25.0 


- 15.4 


09.6 


0.12 


16 A 






1 


| 


20 


0.74 


- 26.7 


- 16.1 


10.6 


0.13 


17 






j 




24 


1.13 


- 39.0 


- 16,1 


22.8 


0,28 


17 A 






1 


i 

j 


24 


1.16 


- 40.5 


- 15.2 


25.3 


0.31 


18 


' 


f 


1 


l 


28 


1.56 


- 44,0 


- 09.2 


34.8 


0.44 


18 A 


1 


1 


28 


1.53 


- 46.0 


- 12.6 


33,4 


0.42 
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TABLE D-I (CON ’T) EXPERIMENTAL DATA 



Dial Speed Voltmeter Voltmeter Voltmeter Measured 
Run Stern Depth Speed (KTS) Zero(V) AVG (V) Diff . (V) Drag(lbs) 



19 2 


1 


20 


0.75 


-44.0 


-28.6 


15.4 


0.19 


19A | 




/ 




20 


0.73 


-45.0 


-29.4 


15.6 


0.20 


20 








24 


1.14 


-46.0 


-20.0 


26.0 


0.32 


20A 








24 


1.16 


-46.5 


-17.9 


28.6 


0.36 


21 








28 


1.61 


-44.5 


-13.2 


31.3 


0.39 


21A 




1 


28 


1.56 


-43.5 


-17.9 


25.6 


0.32 


22 




4 




20 


0.80 


-24.0 


-09.3 


14.7 


0.16 


22A 2 






20 


0.73 


-24.5 


-11.3 


13.2 


0.14 


23 








24 


1.15 


-28.0 


— 09 • 0 


19.0 


0.21 


23A 








24 


1.14 


-26.0 


-07.5 


18.5 


0.21 


24 






f 


28 


1.67 


-25.0 


-04.6 


20.4 


0.23 


24A 




2 


28 


1.60 


-24.5 


-04.0 


20.5 


0.23 


25 




3 


20 


0.73 


-37.0 


-13.2 


23.8 


0.28 


25A 








20 


0.75 


-33.0 


-10.3 


22.7 


0.26 


26 








24 


1.18 


—32 . 0 


-09.4 


22.6 


0.26 


26A 








24 


1.18 


-30.0 


-07.6 


22.4 


0.26 


27 - 








28 


1.53 


-29.0 


-05.5 


23.5 


0.27 


27 A 2 


3 


28 


1.50 


-26 .0 


-03.8 


22.8 


0.26 


28 4 


r 




20 


y0.70 


-40.0 


-08.6 


21.4 


0.25 


28A 








20 


0.70 


-42.5 


-17.5 


25.0 


0.29 


29 








24 


1.00 


-45.0 


-19.2 


25.8 


0.31 


29A 








24 


1.25 


-42.5 


-16.2 


26.3 


0.31 


30 








28 


1.65 


-40.0 


-11.2 


28.8 


0.35 


30A 




r 




28 


1.67 


-40.0 


-11.2 


28.8 


0.35 


31 




2 


20 


0.78 


-35.0 


-17.2 


17.8 


0.20 


31A 




' ' 




20 


0.78 


-34.5 


-15.7 


18.8 


0.21 


32 4 






24 


1.22 


-31.5 


-08.8 


22.7 


0.26 


32A 








24 


1.23 


-32.0 


-10.0 


22.0 


0.25 


33 








28 


1.63 


-32.0 


-07.4 


24.6 


0.29 


33A 




2 


28 


1.58 


-30.0 


-05.2 


24.8 


0.30 


34 




1 


20 


0.77 


-43.0 


-19.7 


23.3 


0.26 


34A 








20 


0.76 


-48.0 


-24.1 


23.9 


0.28 


35 








24 


1.20 


-47.0 


-15.2 


31.8 


0.39 


35A 








24 


1.18 


-40.5 


-12.3 


28.2 


0.34 


36 








28 


1.58 


-41.0 


-05.4 


35.6 


0.45 


36A 4 


1 


28 


1.58 


-44.0 


-07.7 


36.3 


0.45 


37 5 






20 


0.78 


-37.0 


-18.1 


18.9 


0.21 


37A j 




i 




20 


0.80 


-44.0 


-20.9 


23.1 


0.27 


38 








24 


1.21 


-30.5 


-05.0 


25.5 


0.30 


38A 








24 


1.18 


-34.5 


-06.2 


28.3 


0.34 
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TABLE D-I (CONT *D) EXPERIMENTAL DATA 



Dial Speed Voltmeter Voltmeter Voltmeter Measured 
Run Stern Depth Speed (KTS) Zero (V) AVG (V) Diff.(V) Drag(lbs) 



39 






i 


28 


1.67 


-37.5 


-01.1 


36.4 


0.46 


39A 






1 


28 


1.80 


-43,0 


-07.3 


35.7 


0.44 


40 






2 


20 


0.80 


-45.0 


-20.1 


24.9 


0.29 


40A 


5 






20 


0.80 


-46.0 


-20.3 


25.7 


0.30 


41 










24 


1.28 


-46.5 


-15.3 


31.2 


0.38 


4lA 










24 


1.28 


-46.0 


-13.5 


32.5 


0.40 


42 










28 


1.70 


-46.0 


- 8.8 


37.2 


0,47 


42A 






2 


28 


1.65 


-46.0 


- 8.4 


37.6 


0.47 


43 






< 




20 


• 0.75 


-43.5 


-18.1 


25.4 


0.30 


43A 










20 


0.75 


-39.5 


-17.3 


22.? 


0.26 


44 










24 


1,25 


-39.0 


-14.4 


24.6 


0.29 


44A 










24 


1.28 


-38.5 


-13.1 


25.4 


0.30 


45 


5 


f 


' 


i 


28 


1.63 


-37.0 


-07.4 


29,6 


0.35 


45A 


5 


r 




28 


1.60 


-36.0 


-07.0 


29.0 


0.35 



Depth 1 - 18” to Keel (4.5*Dia.) 

Depth 2 - 10” to Keel (2.5*Dia.) 

Depth 3 “ 3.6”to Keel (0.9*Dia.) 

Runs - 1 - 21A Conducted 24 March 1962 
22 - 45A Conducted 25 March 1962 

Water Temp = 51° F 

“I/'* - 1,386 x 10“ 5 ft 2 /sec 

S * 1.939 lb - sec 2 /ft^ 

Dial speed is setting on manual speed control knob. 
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Run 

1 

1A 

2 

2A 

3 

3A 

4 

4A 

5 

5A 

6 

6A 

7 

7A 

8 

8A 

9 

9A 

10 

10 

11A 

11B 

12 

12A 

13 

13A 

14 

14A 

15 

15A 

16 

16A 

17 

17A 

18 

18A 

19 

19A 

20 

20A 

21 

21A 



TABLE D-II - SEPARATION POINT DATA 



Stern Run 

Classification 



3 LD 

k LD 

MD 
MD 
HD 
HD 
LM 
LM 
MM 
MM 
HM 
HM 
LS 
LS 
MS 
MS 

v HS 

3 HS 

1 LS 

LS 
MS 
MS 
HS 
HS 
LM 
LM 
MM 
MM 
HM 
HM 
LD 
LD 
MD 
MD 

v HD 

1 HD 

2 LD 

A LD 

MD 

MD 

HD 

HD 



Separation Body 
Point Angle 



_ ^ H 

5,0 


33 


5.0 


33 


5.0 


33 


5.0 


33 


5.0 


33 


5.0 


33 


4.75 


28 


4.75 


28 


4.75 


28 


4.75 


28 


4.75 


28 


5.0 


33 


4,5 


23 


4.5 


23 


4.5 


23 


4.5 


23 


4.75 


28 


4.5 


23 


0.8 


24 


0.8 


24 


1.25 


42 


1.25 


42 


1.3 


42 


1.35 


42 


0.75 


24 


0.75 


24 


1.0 


32 


1.1 


32 


1.1 


32 


1.1-1. 2 


32 


0.8-09 


24 


0.8 


24 


0.8 


24 


0.8-0. 9 


24 


1.3 


42 


1.3 


42 


2.75 


30 


2.5 


25 


3.0 


36 


3.1 


36 


? 


— 


3.0 


36 
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TABLED-!! (cont ' d) 



Run 


Stern 


Run 

Classification 


Separation 

Point 


Body 

Angle 


22 






LM 


3.o" 


36° 


22A 






LM 


2.6 


25 


23 






MU 


2.6 


30 


23 A 






MM 


2.8 


30 


24 






HM 


? 


— 


24A 






HM 


2.75 


30 


25 






LS 


2.9 


33 


25A 






LS 


2.8 


30 


26 






MS 


3.1 


36 


26A 






MS 


3.2 


42 


27 


> 




HS 


3.3 


42 


27A 


2 


HS 


3.25 


42 


28 


4 


LS 


7.0 


38 


28A 






LS 


7.0 


38 


29 






MS 


7.1 


44 


29A 






MS 


7.1 


44 


30 






HS 


7.3 


49 


30A 






HS 


7.3 


49 


31 






LM 


7.1 


44 


31A 






LM 


7.1 


44 


32 






MM 


7.2 


49 


32A 






MM 


7.2 


49 


33 






HM 


7.2 


49 


33A 






HM 


7.2 


49 


34 






LD 


6.8 


28 


34A 






LD 


6.7 


28 


35 






MD 


6.8 


28 


35A 






MD 


6.9 


-30 


36 






HD 


7.0 


38 


36A 


4 


HD 


? 


— 


37 


5 


LD 


10.5? 


23 


37A 






LD 


? 


- 


38 






MD 


10.5 


23 


38A 






MD 


10.5 


23 


39 






HD 


10.5 


23 


39A 






HD 


10.5 


23 


40 






LM 


11.0 


30 


40A 






LM 


11.0 


30 


41 






MM 


? 


— 


41A 






MM 


? 


— 


42 






HM 


10.9 


30 


42A 






HM 


? 


- 


43 






LS 


? 


- 


43A 






LS 


? 


- 


44 






MS 


11.5 


60 


44A 






MS 


11.5 


60 


45 






HS 


11.5 


60 


45A 


1 




HS 


11.5 


60 
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NOTE : 



TABLE D— II (Cont’d) 



Run classification code - 1st letter designates speed 

L - low 
M - medium 
H - high 

2nd letter designates depth 

D - depth 1 (deep) 

M - depth 2 (medium) 

S = depth 3 (surfaced) 

Separation point is measured in inches from end of cylindrical 
section. Marks on photographs are \ inch apart. 

Body angle is angle between body axis and body surface at 
separation point. 
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APPENDIX E - DISCUSSION OF DRAG FORCE MEASUREMENTS 



The experimental data from Table D-I was reduced to co- 
efficient form using the following procedure. 

(Drag measured) - (Strut drag) = (Total corrected drag) 
where strut drag was calculated by the method that follows. 
Measurement of a strut tare drag was attempted but due to 
either the end effects on the strut because of the lack of 
a dummy body or faulty instrumentation these values were 
of a larger magnitude than the values of measured drag of 
the body plus strut and thus were not valid. An engineer- 
ing value of strut drag composed of profile drag, spray 
drag, wave drag, and interference drag with the body was 
calculated as follows: 

Profile Drag: 

t/„ - .8 - .246 for strut used. 

3.25 

Reynold Number on Strut '* 3 x 10^ 

from figure 126, page 247 of [9] 

Cd - 0.08 

Drag “ q Cd (.8 x SL') - q (.064) SL 
144 144 

where SL - strut length submerged (inches) 

q - 1 ^ V k 2 (1.689) 2 
2 

Spray Drag: 

Froude number on strut ^ 4 x 10 

t =* .8 - .246 

2x 3.25 

from figure 25, page 10-13 of [l0) 
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x * 1 . 625 



Cd x - .24 |'t j 2 - . 058 

Drag - q Cd (X 2 ) - q (.1535) 

X 144 

Wave Drag: 

from figure 24, page 10-13 of (lOj 
Cd t =s 1.4 

Drag “ q Cd + t^ - q (.895) 

144 

Interaction Drag: 

from figure 23, page 8-10 of (lOj 

C ( j t = ,75 x t /c = ,19 

Drag = C d q t 2 = q (.121) 
t lT4 



Total Strut Drag: 

Strut Drag “ 1_ j> V K 



2 



(1.689) 2 F.064SL +1,17 

144 i 



(27) 



The total corrected drag was converted to coefficient 
form by the following relationship: 

Cf “ Drag (28) 

1/2 ? V K (1.689) 2 ( $ ) 
where S = wetted surface of the model. 

The coefficient of frictional resistance, Cj, was neA calcu- 
lated using Schoenherr's relationship that: 

C f = 0.0585 jlog 10 (R h C f )j -2 

« 0,0585 o (29) 

£4343 log e (Rn C f )j 

In order to use the electronic computer to evaluate (29) } 
the Schultz-Grunow formulation 
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Cf - 0.0725 [log 10 Rn -2.0 | 
- _fl , . 07 25 



( 30 ) 



[.4343 log e R n -2. oj 2 
which yields Cf values greater than that of (29) in the range 
of interest in this problem was calculated and used as a first 
value to solve (29) which is implicit in Cf. 

The residual resistance coefficient, C r , was then 
calculated from: 



c t - c f 



(31) 



It was hoped that perhaps some correlation between the 
change of location of separation on the body and change in C r 
for various runs could be found if our drag measurements were 
good. 

An IBM 7090 comuter program, Program I was written to 
solve equations (27) through (31). The details of this program 
are given in Appendix C., Sample outputs of this program are 
included in Table E-I at the end of this appendix. These values 
which are typical of the values obtained on all the bodies are much 
too large for the underwater shapes tested. At submerged depth 
1 the values of Cr and C^ are over ten times as large as would 
be expected by comparison to Series 58 data. 

At first strut drag was suspected but in the surfaced con- 
dition where strut drag does not enter into the calculation the 
residual drag values are much too large also. The bow waves 
generated in the surfaced condition are shown in Figure IX. It 
is not possible for the wave making resistance created in this 
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Besides 



condition to account for such large values of C . 

r 

being too large in magnitude the values are not consistent 
from body to body with the highly streamlined body some- 
times showing a higher value than a less streamlined body 
for similar speed and depth of submergence. Plots of C r 
versus both Reynolds number and speed-length ratio were 
prepared but these were completely inconsistent and unsatis- 
factory when compared to similar plots of similar bodies in 
the literature. 

At the time of this writing the cause of the excess 
drag is unknown. It is speculated that the dynamic response 
of the drag measuring system is causing a false reading of 
drag due to some phenomena that does not appear during sta- 
tic calibration. The calibration curves used in this work 
are shown in Figure E-I. It is also of interest that a 6 - 
nother group using another submerged model experienced the 
same difficulty, i.e., too much drag. 

If it is desired to make further drag measurements 
on the model it is recommended that the drag force dyna- 
mometer be mounted at the model rather than in the strut. 
This will eliminate the problem relating to strut drag with 
the exception of the interference drag between the strut and 
body. It is also recommended that the dynamic response of 
the entire drag measuring system used during the tests be 
given a thorough investigation. The authors realize that 



87 



the new towing tank and dynamometer equipment was put to 
use on an emergency basis, the construction of such appara 
tus being completed, and before the carriage and instrumen 
tation were properly calibrated, accepted, and put into 
regular use. Upon completion of all of these steps it 
should be feasible to carry out successful drag measure- 
ments on submerged bodies at the new M.I.T. Towing Tank. 
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Figure E-r I- Dynamometer Calibration 




Volts (£>.c.) 



TABLE E-I - SAMPLES OF COMPUTER OUTPUT, PROGRAM I 



BCCY f\iUN , fcfcR= l 



UEPTh= 1 



V 6 L • (KT. 


RLYNOLUS NO. 


SPE EO/L ENGTH 


TOT. DRAG 


TOT.URA 


0.74 


2. 3726 


0 5 


C • 4 6 


0. L20 


C.098 


C • 74 


2.372L 


05 


C . 46 


0. L30 


0. 108 


i. 13 


3.6226 


05 


C . 70 


0.280 


0.229 


1. L6 


3 . 7 1Q6 


05 


0.72 


0.3L0 


0.257 


L • 56 


5. CGOt 


05 


C . 96 


0.440 


0.343 


L • 5 3 


4.9046 


05 


C • 94 


0.420 


0.327 



DEPTh=2 



VEL • ( KT . 


REYNOLDS NO . 


SPEED/LENGTH 


TOT. DRAG 


TOT. DRAG 


0.75 


2.4C4E 


05 


C . 46 


C. L60 


0.L43 


0.75 


2 . 404t 


05 


0.46 


C. L50 


0. 133 


L. L 9 


3.8L46 


0 5 


C . 7 3 


0.280 


0.238 


L. LL 


3 • 5 58t 


05 


0.68 


0.280 


0.243 


L.6L 


5. L60E 


05 


C . 99 


0. 38C 


0.303 


1.6L 


5. L 606 


05 


0.99 


0.390 


0.313 



06 P Tl-= 3 



VEL. (KT. 


REYNOLDS NO. 


SPE EO/LENGTH 


TOT. DRAG 


TOT. DRAG 


0.77 


2 • 468t 05 


0.47 


0.280 


0.280 


0.77 


2.4686 05 


0.47 


0.310 


0.310 


1.16 


3.718E 05 


C . 72 


0.280 


0.280 


1.21 


3.8786 05 


0.75 


0.260 


0.260 


1.58 


5 • C 64 6 05 


C.97 


0.330 


0.330 


1.58 


5 • C64t 05 


0.97 


0.320 


0.320 




CT 


CF 


CR 






2.541E-02 


5 • 909E-0 3 


1.950E-02 






2.799E-02 


5 • 909E-03 


2.208E-02 




DEPTH 


= 1 2.543E-02 


5.400E-03 


2 .0036-02 






2 • 700E-02 


5.371E-03 


2. 163E-02 






1.998E-02 


5.053E-03 


l .49 36-02 






1.978E-02 


5.073E-03 


1.471E-02 





DEPTfl * 2 



CT 

3.604E-02 

3.353E-02 

2.377E-02 

2*7956-02 

1.653E-02 

1.707E-02 



CF 

5.892E-03 
5 • 892E-03 
5.342E-03 
5.421E-C3 
5.021E-03 
5.021E-03 



2.763E-02 
1.842E-02 
2.253E-0 2 
1.L5LE-02 
1.205E-0 2 



DEPTH 



CT 

7. 797E-02 
8.632E-02 
3.436E-02 
2.932E-02 
2.182E-02 
2. LL6E-02 



CF 

5. 858E-03 
5.8586-03 
5.37LE-03 
5.324E-03 
5.040E-C3 
5.040E-03 



CR 

7.2LLE-02 
8.047E-02 
2 • 898E-02 
2.400E-02 
L.678E-02 
L.6L2E-02 



90 



TABLE E-I (CON’T) 



BODY NUNBER= 


3 








D E P T F = 1 










VEL. ( KT . 


RLYNLLDS NO. 


SPEED/LENGTH 


TOT. DRAG 


TOT. DRAG (CORK) 


0. 80 


2.856E 05 


0.47 


0.210 


0.185 


0.77 


2.749E 05 


C • 4 5 


0.2 10 


0. 186 


1. 12 


3.999E 05 


0.65 


0.330 


0.280 


1.1 2 


3.999B C5 


C . 6 5 


0.340 


0.290 


UAb 


5.2 84t 05 


C . 86 


0. 380 


0.293 


1.48 


5.284b 05 


C • 86 


0.400 


0.313 


DEPTH=2 










VEL. ( KT. 


REYNOLDS NO. 


SPEED/LENGTH 


TUT. DRAG 


TOT. DRAG (CORR) 


0.75 


2.678E 05 


C • 44 


0.200 


0.183 


0.74 


2.642b 05 


0.43 


0.210 


0.194 


1.09 


3.892b 05 


0.64 


0.290 


0.255 


1.C9 


3.892E 05 


C . 64 


0.330 


0.295 


1.5C 


5 • 356t 05 


0.88 


0.310 


0.243 


1.51 


5.392E 05 


C . 88 


0.310 


0.242 


DEPTHS 










VEL. (KT. 


REYNOLDS NO. 


SPEED/LENGTH 


TOT. DRAG 


TOT. DRAG (CCRR) 


0.77 


2 • 749t 05 


0.45 


0. 190 


0. 190 


0.79 


2.821t 05 


C . 46 


0. 190 


0.190 


1 « 1 & 


4-.-14 2E 05 


0.68 


0.270 


0.270 


1.17 


4.178t 05 


C • 68 


0.250 


0.250 


1.62 


5.784b 05 


0.95 


0.330 


0.330 


1.58 


5.641E 05 


C.92 


0. 300 


0.300 




CT 


CF 


Ck 






3.771E-02 


5.678E-03 


3.203E-02 






4. U1E-02 


5.724E-03 


3.539E-02 




-DEPTH 


= 1 2 • 920E-02 


5 • 290E-03 


2.391E-02 






3 • 024E-02 


5.290E-03 


2.495E-02 






1. 749E-02 


4.997E-03 


1.249E-0 2 






1.868E-02 


4.997E-03 


1.369E-02 





CT 


CF 


CR 


4.258E-02 


5.756E-03 


3.682E-02 


4.623E-02 


5.773E-03 


4.046E-02 


DEPTH - 2 2.801E-02 


5. 320E-03 


2.269E-02 


3.241E-02 


5.320E-C3 


2.709E-02 


1.411E-02 


4.983E-03 


9. 126E-03 


1.307E-O2 


4.976E-03 


8.894E-03 





CT 


CF 


CR 




5.058E-02 


5. 724E-03 


4.485E-02 




4.805E-02 


5.693E-02 


4.235E-02 


DEPTH =3 


3. 167E-02 


5.252E-03 


2.642E-02 




2.B82E-02 


5.243E-03 


2.358E-02 




1.985E-02 


4.906E-03 


1.494E-02 




1.897E-02 


4.931E-03 


1.403E-02 
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TABLE E-I (CON*T.) 



BODY NUKBER=5 



OEPTF= 1 



VEL. (KT. 


REYNOLDS NO. 


SPEED/LENGTH 


TOT. DRAG 


TOT. DRAG 


0.78 


3 . 24 IE 05 


0.42 


0.210 


0.186 


0.80 


3.324E 05 


0.43 


0.270 


0.245 


1.21 


5.028E 05 


0.66 


0.300 


0.242 


1.18 


4.903E 05 


0.64 


0.340 


0.285 


1.67 


6 • 94Gt 05 


0.90 


0.460 


0.349 


1.8C 


7.480E 05 


C . 97 


0.440 


0.311 



OEPTH=2 



VEL. (KT. 


REYNOLDS NO. 


SPEED/LENGTH 


TOT. DRAG 


TOT. DRAG 


0.80 


3.324E 05 


0.43 


0.290 


0.271 


0.80 


3.324b U 5 


0.43 


0.300 


0.281 


1.28 


5.319F 05 


0.69 


0.380 


0.331 


1.2b 


5.319E 05 


C . 69 


0.400 


0.351 


1.70 


7.064b 05 


0.92 


0.470 


0.384 


1.65 


6.857b 05 


C • 89 


0.470 


0.389 



DEPTH=3 



VEL. (KT. 


REYNOLDS NO. 


SPEED/LENGTH 


TOT. DRAG 


TOT. DRAG 


0.75 


3. 1 17t 05 


0.41 


0.300 


0.300 


0.75 


3.117E 05 


0.41 


0.260 


0.260 


1.25 


5.194E 05 


C • 68 


0.290 


0.290 


1.2b 


5 . 3 19t 05 


0.69 


C.300 


0. 3C0 


1.63 


6. 773t 05 


C . 88 


0.350 


0.350 


1.60 


6.649t 05 


0.87 


0.350 


0.350 




CT 


CF 


CR 






3. 540E-02 


5 . 528E-03 


2.987E-02 






4.4296-02 


5.498E-03 


3 . 879E-02 




DEPTH 


= 1 1.915E-02 


5.047E-03 


1.410E-02 






2. 3706-02 


5.073E-03 


1.863E-02 






1.4526-02 


4.731E-03 


9.785E-03 






1.1 146-02 


4.662E-03 


6.478E-03 





-DEPTH 



CT 

4-905E-02 
5.086E-02 
2.342E-02 
2.483E-02 
1.539E-02 
1 • 655E-02 



CF 

5.498E-03 
5.498E-03 
4 • 990E— 03 
4 • 990E-03 
4.714E-03 
4.742E-03 



CR 

4.355E-02 
4.537E-02 
1 • 843E-02 
1.984E-02 
1.067E-02 
1. 181E-02 



DEPTH 



CT 

7.374E-02 
6.391E-02 
2.566E-02 
2 • 5 32E-02 
1.821E-02 
1 • 890E—02 



CF 

5 • 574E-03 
5 • 574 E-03 
5.014E-03 
4 • 990E-03 
4.754E-03 
4.771E-0; 



CR 

6.817E-02 
5.834E-02 
2.065E-02 
2.033E-02 
1 • 346E-02 
1.413E-02 
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APPENDIX F - SAMPLE CALCULATIONS 



Calculation of Hydrostatic Pressure. 

The hydrostatic pressure to the keel of the bodies at 
the three depths of submergence is desired for use in com- 
puter Program III to compute the gage pressure on the body 
from the following relationship; 

Pgage - /P\Q + Ps (144) 



(i)> 



where 



4 

P - dynamic pressure (ps£) 

% - 1/2 <^LC 2 (psf ) 

Ps “ hydrostatic pressure (psi) 



Ps 



59 h 

1725 



where <? - 1.935 lb.sec^ at SI 0 ! 1 

9 - 32.174 ft /sec 2 
h - depth to keel (inches) 

Using the foregoing relation the following results were 
obtained: 



Depth 



h (inches) 



Ps (psi) 



1 

2 

3 



18 
10 -• 
3.6 



.6484 

.3602 

.1297 
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Calculation of Point of Separation from Stratford Criterion 

A typical calculation of the location of the separation 
point according to Stratford’s criterion is given below for body 
#3 at velocity of 0.75 knots. 

From the form of (21) separation can occur only when 
Cp , as plotted in Figure F-II, is positive since negative 
values of Cp will result in imaginary numbers. 

At X - 34.8” X - 34.8 -0.99 
L 35.15 

from Figure F-I; U - 0.75 

Uoo 

from Table (o-I; Reynolds number - 2.65 x 10® at X-34.8” 

thus its = .75 (2.65 x 10®) - 1.99 x 10® 

from Figure F-II; Cp = 0.05 and dCp - 0.50 

dx 

n - log 10 Rs - 5.30 from (22) 
o 

since, dp £ o; /3 - .73 from (24) 
dx^ 

thus the left hand side of (21) becomes 

5.30-2 1 

(2x.05)““? (34.8 x .50) 2 = .625 

while the right hand side of (21) becomes 

• „5 a “1/ in 

(1.06) (.73) (1.99 x 10 xl0 D ) - .912. Thus the 

criterion is not satisfied at X “• 34.8" since the RHS is much 

larger than LHS. Next try 

X - 34 . 9 X - .995 
L 



94 



0.96 



from Figure F-II; Cp - 0.10 and dC 






thus the LHS of (21) becomes 




1.46 which now makes the 



RHS smaller than LHS so the point of separation lies between 
X - 34.8" and X - 34.9." 

Calculations similar to the above were used to cal- 
culate the separation points on the other bodies at the va- 
rious speeds. During these calculations it was noted that the 
criteria was insensitive to speed changes and that the same 
point of separation occurred for all values of speed. The 
points of separation according to this criterion are marked 
by small circles on Figure F-II. 

Calculation of Point of Separation from Buries Criterion 

A typical calculation of the location of the separation 
point according to Buri's criterion is given below for body #3 
at a velocity of 1.12 knots* 

From the form of equation (26) separation will occur 
only when <1 / U\ as determined from Figure F-I is negative. 



Uoo 




at X - .906; 

E 




—.115 and U » 1.103; from Figure F-I 
Uoo 



and 0 - .091; from Table G-I 
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thus the left hand side of (26) becomes 



-(,115)(1.103)(. 906X091) - -.0105 > — .006 

Thus the point of separation lies slightly forward of X - .906 

L 

but in the region of where 

d/_U_\ is still negative. Separation point is at 

uooy 



the point can be calculated within the accuracy of the plot in 
Figure F-I . 

Calculations similar to the above were used to calculate 
the separation points on other bodies at various speeds. During 
these calculations it was noticed that the criterion was insensitive 
to speed changes and that the same point of separation occurred 
for all values of speed. The points of separation according to 
this criterion are marked by small circles on Figure F-I. 



Figures F-I and F-II on the next two pages were plotted 

from values obtained from the outputs of Programs II and III 

as aids to evaluate (21) and (26). 




X « .900 or X = 31.7 inches which is as close as 

E 



NOTE: 
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FIGURE F-I - PLOT OF U vs X 

TJoo L 




L 
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FIGURE F-II - Plot- of ys X (distance from body nose) 
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APPENDIX G - SAMPLES OF COMPUTER OUTPUT. PROGRAMS Hand IP 



TABLE G-I - SAMPLES OF COMPUTER OUTPUT, PROGRAM XI 



BODY NUMBEK=3 
BOGY DEPTHS l 
VELOCITY(KTS)=i. I? 



X ( INCH) 


RADUI S 


CORK. RADIUS 


DELTA 


DELTA-STAR 


THE TA 


REYNOLDS NO. 


Y-NORM/ 


0,19 


0.40 


0.401 


.0154 


.0019 


.0015 


2. 161E 


03 


0. 22 825 E- 01 


0.97 


0.94 


0.946 


.0564 


.0068 


.0056 


1. 103E 


04 


0. 53822E-0 1 


' 2.36 


1.37 


1.383 


. 1141 


.0135 


". 0114 


2 .68 4 E 


04 


0. 78683E-0 1 


4.30 


1.67 


1.691 


. 183? 


.0214 


.0184 


4.891E 


04 


0. 96 202 E- 0 1 


6.72 


1.87 


1.900 


.2599 


.0298 


.0263 


7.643E 


04 


0. 10806E-00 


9.53 


1.97 


2.003 


. 3408 


.0383 


.0348 


1.084E 


05 


0. 11424E-0Q 


12.62 


2.00 


2.047 


.4225 


.0465 


.0437 


1.435E 


05 


0. 1 1 64 l E-00 


15.90 


2.00 


2.054 


.5032 


.0542 


.0526 


1.808E 


05 


0. U685E-00 


19.24 


2.00 


2.06L 


.5806 


.0613 


.0614 


2. 188E 


05 


0. 1 17 2 5 E - 0 0 


22.52 


2.00 


2.066 


.6528 


.0677 


.0697 


2.561E 


05 


0. U762E-00 


2 5.6 2. 


2.00 • 


2.073 


.7181 


.0734 


. 0773 


2.914E 


05 


0. 11794E-00 


26.43 


2.00 


2.078 


. 775 1 


.078? 


. 084 1 


3.234E 


05 


0. 1 1821E-00 


30.85 


1.94 


2.C21 


.8206 


.0614 


.0899 


3.509E 


05 


oTi 1497E-00 


32.79 


1.6? 


1.697 


• 8_42 7 


.080? 


. 0947 


3. 729E 


05 


0.96544E-01 


34.17 


1 . 15 


1.217 


.*8326' 


.07 34 


.0985 


3.886E 


05 


0 • 692 l 3 E-0 1 


34.96 


U . 57 


0.605 


.7529 


.0560 


. 1022 


3.976E 


05 


0. 34408E-01 



BODY NUMBER=3 
BOGY L.EPThs? 



VELCCITY(K1 S ) = 1 • 5 1 



X [ INCH) 


KADUIS 


CORK. RADIUS 


DELTA 


DELTA-STAR 


THETA 


REYNOLDS NO. 


Y-NORM/ 


0. 19 


0.40 


0.401 


.0145 


.0018 


.0014 


2.914E 


03 


0. 22 82 1 E-0 l 


0.97 


0.94 


0.946 


.0531 


.0064 


.0053 


1.487E 


04 


0. 5 3R03E-0 1 


2.36 


1.37 


1.383 


. 1076 


.0128 


.0107 


3.619E 


04 


0. 78642E-0 l 


4.30” 


1.67 


1.690 


.1728 


.0202 


.01 73 


6 . 5 9 4 E 


04 


0 . 96 1 38 E-0 1 


6.72 


1.87 


1.898 


.2453 


.0282 


.0248 


1.030E 


05 


0. 10797E-00 


9.53 


1.97 


2.006 


.3218 


.0364 


.0328 


1.461E 


05 


0. 1 1413E-00 


12.62 


2.00 


2.044 


.3991 


.0442 


.0411 


1.935E 


05 


0. 1 1628E-00 


15.90 


2.00 


2.052 


.4756 


.0516 


.0495 


2 • 4 3 8 E 


05 


0. 1 1670E-00 


19.24 


2.00 


2.058 


.5490 


.0586 


. 0578 


2.950L 


05 


0. 1 1709E-00 


22.52 


2.00 


2.065 


.6175 


.0646 


.0656 


3.453E 


C 3 


0. 1 1 7446-00 


23.62 


2.00 


2.070 


.6795 


.0701 


.0728 


3. 929E 


05 


0. 1 1 775 E-00 


28.43 


2.00 


2.075 


.7333 


.0747 


.0792 


4. 360E 


05 


0. 1 1 80 1 E-00 


30.83 


1.94 


2.018 


.7771 


.0779 


.Of 46 


4.731E 


05 


0. 11477E-00 


32.79 


1.62 


1.694 


.7987 


.0769 


.0891 


5.028E 


05 


0.96362E-01 


34.17 


1.15 


1.214 


.7903 


.0706 


.0927 


5.240E 


05 


0. 69068E-0 l 


34l96 


0.57 


0.604 


.7173 


.054? 


.096 1 * 


5.361E 


05 


0. 34 34 3E-0 1 



99 










I 



I 

1 



I 



I 

I 













ROC Y NUMB EK = l 



80CY DEPTHS? 
VEL0CITY(KTS)=0.75 



X< INCH) 


RAUUIS 


CORR. RADIUS 


DELTA 


DELTA-STAR 


THETa 


REYNOLDS NO. 


Y-NORM/ 


0.17 


0.3 7 


0. 371 


.0152 


.0019 


.0015 


1.295E 


0 3" 


0. 23524E-0 l 


0.87 


0.88 


0.R8o 


.0559 


.0067 


.0055 


6.626E 


03 


0. 56149E-01 


2. 12 


i .30 


1.309 


.1134 


.0134 


.0113 


1.615E 


04 


0. 82974E-0 1 


.3.8*3 


L .60 


1.621 


. 1015 


.0211 


.0182 


2.932E 


04 


0. 10272E-00 


6.02 


L .81 


1.839 


.2576 


.0295 


.0261 


4. 585E 


04 


0. 1 1656E-00 


R. 53 


L .96 


L. 998 


. 3379 


.0380 


.0346 


6.497E 


04 


0. 12661 E- 00 


11.31 


2.00 


2 . 0 ■ ♦ 6 


.4195 


.0462 


.0433 


8.614E 


04 


0. 12967E-00 


14. 2S 


2.00 


2.054 


.4997 


.0539 


. 0522 


1.085E 


05 


0. 13016E-00 


17.24 


2.00 


2.061 


.5765 


.0610 


.0609 


1.313E 


05 


0. 13061E-00 


20.19 


2.00 


2.067 


.6484 


.0674* 


.0692 


l. 538E 


05 


0. 13 101E-00 


22.96 


2.00 


2.073 


.7131 


.0729 


.0767 


1.749E 


05 


0. 13137E-00 


25.48 


2.00 


2.078 


.7699 


.0777 


. 0835 


1.941E 


05 


0. 1 3 167E-00 


27.65 


2.00 


2.0R2 


.8173 


.0816 


.0892 


2. 1 06E 


05 


0. 13191E-00 


29.38 


2.00 


2.C85 


.8543 


.0846 


.0937 


2.238E 


05 


0. 1 32 10 E -0 0 


30.63 


1.63 


1.69 7 


.8628 


.0818 


. 0972 


2.333E 


05 


0. 10754E-00 


31.33 


0.78 


0.805 


.7922 


.0637 


. 1005 


2. 386E 


05 


0. 50993E-0 l 



BODY NUMbEK=3 
BODY L) E P T H ■= 2 



VELOC I T Y < KTS ) =0 • 75 



X{ INCH) 


RAUUIS 


CCRR . RADIUS 


delta 


DELTA-STAR 


THETA 


REYNOLDS NO. 


Y-NORM/ 


0. 19 


0.40 


0.401 


.0167 


.0020 


.0016 


1.447E 


03 


0. 22831E-01 


0.97 


0.94 


0.947 


.0610 


.0073 


. 0060 


7.388E 


03 


0. 53850E-0 l 


2.36 


1.37 


l . 384 


. 1234 


.014o 


.0123 


1 • 797 E 


04 


0.7874 1 E-0 1 


4 . 3 U 


1.67 


1.693 


. 1980 


.0230 


.0199 


3.275E 


04 


0 . 96293E-0 l 


6.72 


1.87 


1.902 


.2808 


.0320 


.0285 


5. 118E 


04 


0. 108 19E-00 


9.53 


l .97 


2.011 


. 3679 


.0410 


.0378 


7.258E 


04 


0. 11439E-00 


12 . 6 ? 


2.00 


2.050 


.4559 


.0497 


.0474 


9.612G 


04 


0. 1 1659E-00 


15.90 


2.00 


2.058 


.5426 


.0579 


.0570 


1.211E 


05 


0. 11706E-00 


19.24 


2. CO 


2.065 


.6256 


.0654* 


~. 0665 


~i;465E 


05 


o. “ii7 * 5 bf - o'er 


22.52 


2.00 


2.072 


.7030 


.0721 


. 0756 


1.715E 


05 


0. 11 787E-00 


25.62 


2.00 


2 . C73 


.7730 


.0780 


.0839 


1.951E 


05 


0. 1 1 820E- 00 


2R.43 


2.00 


2.08 3 


.8340 


.0830 


.0912 


2. 165E 


05 


0. 1 1849E-00 


30.85 


1.94 


2.026 


.8820 


.0863 


.0975 


2 • 350 E 


05 


0. 11524E-00 


32.79 


1.6? 


1.702 


.9051 


.0848 


. 1027 


2.497E 


05 


0.96797E-01 


34.17 


1.15 


1.220 


.8923 


'.0772 


. 1069 


2.603E 


05 




34.96 


0.57 


0.606 


.8028 


.0586 


.1110 


2. 663E 


05 


0.34498E-01 



80DY NUMBERS 
80DY DE PTh = 3 



VEL0CITY(KTS)=0.75 



X( INCH) 


RAOUIS 


CORR. RADIUS 


DELTA 


DELTA-STAR 


THETA 


REYNOLDS NO. 


Y-NORM/ 


0.22 


0.42 


0.422 


.0187 


.0023 


.0018 


1.676E 


03 


0. 20604E-01 


1.13 


1.00 


1.008 


.0639 


.0082 


.0068 


8.606E 


03 


0. 492 50E-0 l 


2.75 


1.42 


1.436 


.1392 


.0163 


.0139 


2.094E 


04 


0. 70 191 E- 01 


5.00 


1.72 


1.746 


.2228 


*.02 57 


7 0*2 2 5 


808E 


04 


~0 .8531 6 E-0 1 


7.81 


1.91 


1.946 


.3155 


.0356 


.032? 


5.948E 


04 


0. 95092E-01 


11.08 


2.00 


2.046 


.4130 


.0456 


.0426 


8.439E 


04 


0. 99979E-01 


14.69 


2.00 


2.055 


.5113 


.0550 


.0535 


l. 119E 


05 


0. IC044E-00 


18.50 


2.00 


2.064 


.6076 


.0638 


.0645 


1.409E 


05 


0. 10087E- 00 


22.39 


2.00 


2.072 


. 7000 


.0718 


.0752 


1.705E 


05 


0. 10 1 26E-00 


26. 20 


2.00 


2.079 


.7857 


.0790 


.0854 


1.995E 


05 


0. 10161E-00 


29.81 


1.98 


2.065 


.8626 


.0851 


.0948 


2.270E 


05 


0. 10 093E-00 


33.08 


1.85 


1.938 


.9239 


.0887 


. 1033 


2.519E 


05 


0. 94742F-01 


35.89 


1.57 


1.657 


.9614 


.0882 


. 1106 


2. 7 34E 


05 


0. 8 l OOOE-O l 


38.14 


1.18 


1.260 


.9655 


.0824 


. 1169 


2. 905E 


05 


0. 61602E-01 


39. 75 


0. 73 


0.795 


.9154 


.0692 


. 1225 


3.027E 


05 


0. 38848E-0 1 


40.67 


0.28 


0.316 


.7423 


.0443 


. 1297 


3.098E 


05 


0.1T4TDE-01 
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TABLE G-II - SAMPLES OF COMPUTER OUTPUT, PROGRAM III 



tjonv NUMB EK= 1 



CEPTH=2 



SPEtU=0.75 



X ( I ) 

-0 . 98940 E 00 

— 0 • 944586 00 
-0.86563E 00 

— 0 • 75540E 00 
-0.61788E 00 
-0.45802E-G0 
-0 • 28 160E-00 
-0.95013E-01 

0.95013E-01 
0*28 160E-00 
0 .458C2E-00 
0 • 6 1 78 8E 00 
0.75540E 00 
0 • 8b 56 3E 00 
0.94458E 00 
0.98940E 00 



YU) 

0.23524E-C1 
0 . 56 14 9 E- C 1 
0.82974E-C1 
C. 1C 2 /2E-C0 
0. 1 1686E-C0 
0. 12661E-C0 
0. 12967E-C0 
0. 130166-C0 
0 • 130blE-C0 
0. 13101E-C0 
0. 131 37E-C0 
0. 131o7E-C0 
0. 131916-CO 
0. 13210E-C0 
0. 10754E-C0 
0.50993E-C1 



U ( I ) 

0.66402E GO 
0.986G8E OG 
0- 74 39 3 E 00 
0. 10463E 01 
0.1035CE 01 
0.10413E 01 
0.10222E 01 
0.10106E 01 
0.1009CE 01 
0. 100996 01 
0.10136E 01 
0* 10200E 01 
0.10552E 01 
- 0. 12431E 01 
0. 11424E 01 
0 • 55872E 00 



P ( 1 ) 

0.55907E 00 
0. 276466-01 
0.43910E-00 
-0. 94643E-01 
-0. 72811E-01 
-0. 84395E-01 
-0 • 44932E-0 1 
-0, 21251E-01 
-0. 17995E-01 
-0, 19858E-01 
-0. 27446E-01 
— 0 • 404 7 2E— 0 1 
-0. U343E-00 
— 0 * 54533E 00 
-0. 30503E-00 
0.68783E 00 



PCC( I ) 
0.52739E 02 
0 . 5 19 1 2E 02 
0.52552E 02 
0.51722E 02 
0.51 756E 02 
0.51738E 02 
0.51799E 02 
0.51836E 02 
0.51341E 02 
0.51838E 02 
0.51826E 02 
0 • 5 1 806E 02 
0.51692E 02 
0.51020E 02 
0.51394E 02 
0.52939E 02 



V EL ( I ) / 
0.49802E-00 
0.73956E 00 
0 • 56 1 70E 00 
0 • 78469E 00 
0.77682E 00 
0.781G1E 00 
0.76666E 00 
0* 757936 00 
0.75672E 00 
0*757416 00 
0.76022E 00 
0.76503E 00 
0.79140E 00 
0 . 93233E 00 
0* 85678E 00 
0*4 1904E— 00 



HOOY NUI“!t<E< = 3 








L’E P TH = 2 








SPEED=0. 75 








X ( I ) 


Y( I ) 


U ( I ) 




-0.98940E 00 


0. 22831 E-Cl 


0.701 15E 


00 


-0.94458E 00 


0. }383CE-C1 


0. 100426 


01 


-0.86563E 00 


C. 7&741E-C1 


0. 10542E 


01 


-0.75540E 00 


0.9^293E-Cl 


0. 10473E 


01 


-0.61788E 00 


0. 1 0 8 l 9E-C0 


0. 10410E 


01 


-O.45802E-C0 


C. 1 1439E-C0 


0. 1C281E 


0 1 


-0 . 20 160E-00 


0. 1 1659E-CC 


0. 101666 


01 


-0.95013E-0L 


0 . 1 1 706E-C0 


0. 100886 


01 


0.93'013E-01 


0. 1 1748E-C0 


0. 10081E 


01 


0.28 160E-00 


C. 1 17e7E-C0 


0. 10123E 


01 


0.45602E-00 


C. 1 L72CE-CO 


0. 10082E 


01 


0.61788E 00 


C.11849E-C0 


0. 10435E 


01 


0. 75540E 00 


0. 1 1524E-C0 


0. 1 1069E 


01 


0.86563E 00 


C . 96797E-C 1 


0. 1 0 7 3 7 C 


01 


0.94458E 00 


0.69414F-C1 


0. 10298E 


01 


0 . 9 894 OE 00 


C . 34498E-C 1 


0. 7 7644 E 


CO 



p ( I ) 


PCCt I 


) 


V E L ( I )/ 


0.50839E 00 


0.52660E 


02 


0.52586E 


00 


-0. 85C96E-02 


0. 51856E 


02 


0. 75318E 


00 


-0. 11125E-00 


0. 51696E 


02 


0. 79062E 


00 


-0. 96933E-01 


0. 51718E 


02 


0 • 7 8 55 IE 


00 


-0. 83723E-01 


0.51739E 


02 


0. 78077E 


00 


-0. 56907E-C1 


0.51 780E 


02 


0.77104E 


00 


-0. 33406E-01 


0. 51817E 


02 


0. 76242E 


00 


-0. 17689C-01 


0 • 5 1 84 IE 


02 


0. 75660E 


00 


-0. 16231F-01 


0 . 5 1 844 E 


02 


0. 75606E 


00 


-0. 24782E-01 


0.51830E 


02 


0.75924E 


00 


-0. 16448E — 01 


0 . 5 1 84 3 E 


02 


0. 75614E 


00 


-0. 88884E-01 


0. 51731E 


02 


0. 78262E 


00 


-0.22517E-00 


0.51518E 


02 


0.83016E 


00 


-0. 15292E-00 


0.51631E 


02 


0.805316 


00 


-C. 60430E-01 


0.51775E 


02 


0. 77233E 


00 


0. 397 14E-00 


0.52487E 


02 


0. 58233E 


00 
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BODY NUMB £K = 3 



OEPTH=l 



SPEED=1. 12 



X { I ) 




Y { 1 ) 


U { 1 ) 


-0.98980E 


CO 


0 • 22 825E-0 1 


0.70131E 


-0.9885RE 


00 


0.53822 E-Ol 


0. 10083E 


-0 . 86563E 


GO 


0.76663E — Cl 


0. 1 0 582 E 


-0.75580E 


00 


0.9 6202E — Cl 


0. 10873E 


-0.61788E 


00 


0. 10 8 06 E — CO 


0. 10809E 


-0.85802E- 


'00 


0. 1 L828E-C0 


0. 1 02 HOfc 


-0.28160E- 


GO 


0. 1 168 1E-C0 


0. 1C165E 


-0.9501 3E- 


■01 


0. 11685E-00 


0. 10087 E 


0.98013E- 


•01 


0. 11728E-00 


0. 1C079E 


0.26160E- 


■00 


0. 1 1 762E-C0 


0. 10096E 


0.85802E- 


•00 


0. 1 1 798E-C0 


0. 10189E 


0.61788E 


00 


0. 11321E — CO 


0. 10399E 


0 • 7558 OE 


00 


0. 1 1897E-CC 


0. 1 1068E 


0.86563E 


00 


0. 96588E-01 


0. 107 32E 


0 • 988 58E 


00 


0.692 13E-C 1 


0. 10289E 


0.98980E 


00 


0. 38808E-C1 


0. 77605E 



BODY MJHB Ek = 3 
DEPTH=2 



SPEED= 1.51 



X( i ) 




Y { I ) 


U( 1 ) 


-0.98980E 


00 


0.22821E-CI 


0.70183E 


-0.98858E 


00 


0.53803E-01 


0. 10083E 


-0.86563E 


00 


0. 78682E-01 


0. 10582E 


-0.75580E 


00 


0.96138E-01 


0. 10873E 


-0.61788E 


00 


0. 1C797E-C0 


0. 10809E 


-0 • 85802 E 


-00 


0. 11813E-00 


0. 1C2B0E 


-0.28160E 


-00 


0. 11628E-00 


0. 10165E 


-0.95013E 


-01 


0.11 670E-G0 


0. 10086E 


0.95013E- 


-01 


0.117 09F-C0 


0. 10079E 


0.28160E 


-00 


0. 11788E-C0 


0. 10095E 


0 . 8 5 80 2 E 


-00 


0. 11775E-C0 


0. 10188E 


0.61788E 


00 


0. 1 1801E-C0 


0. 10398H 


0.75580E 


00 


0. 11877E-C0 


0. 1 1062E 


0.86563E 


00 


0 • 96 362E-C 1 


0. 10729E 


0.98858E 


00 


0 • 69068E-C 1 


0. 10288E 


0.98980E 


00 


0.38383E-01 


0. 77580E 



BODY NUMBER 3 5 
DEPTH- 3 



SP EED=0 • 7 5 



X ( I ) 




Y ( I ) 


U( I ) 


-0.98980E 


GO 


0. 206U8E-01 


0.70371E 


-0.98858E 


00 


0.89250E-C 1 


0. 10285F 


— 0. 8t>563E 


00 


0.70191E-01 


0. 10869E 


-0.75580E 


00 


0.85316E-C 1 


0. 10835F 


-0.61788E 


00 


0.95092E-C 1 


0. 10351E 


-0.85802E 


-00 


0.99979E-C1 


0 • 1 02 67 E 


-0.28160E 


-00 


0. 10088E-C0 


0. 10093E 


-0.95013E 


-01 


0. 1C087E-C0 


0. 10071E 


0.95013E 


-0 1 


0. 10126E-C0 


0. 10076E 


0.28160E 


-00 


0. 10161E-C0 
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